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ABSTRACT 


A  general  discussion  of  free  electron  lasers  with  variable 
parameter  wlgglers  Is  presented  with  a  view  towards  their  potential  for  the 
production  of  high  power  optical  radiation  at  reasonable  efficiency.  The 
theoretical  analysis  Is  based  upon  a  one  dimensional  relativistic 
Hamiltonian  formulation  and  Is  developed  in  a  manner  to  take  advantage  of 
the  analogy  between  the  free  electron  laser  process  and  radio  frequency 
accelerators.  Three  promising  operational  modes  are  Identified  and 
analyzed.  The  first  may  be  th'ought  of  as  an  electron  decelerator  and  is 
thought  to  have  the  most  promise  for  single  pass  devices.  Both  oscillator 
and  amplifier  configurations  are  studied.  The  second  Is  based  upon 
adiabatic  trapping  and  detrapping.  Intended  to  reduce  the  spread  In 
electron  energy  typically  induced  by  the  FEL  process.  The  third  Is  based 
upon  the  method  of  phase  area  displacement.  It  has  the  advantage  of  wide 
gain  bandwidth  and  small  induced  energy  spread,  and  Is  thought  to  have  the 
most  promise  for  storage  ring  applications.  Generally  speaking,  it  Is 
found  that  high  peak  power  Is  Intrinsic  to  these  modes  of  operation. 
Potential  problems  from  parasitic  oscillations  analogous  to  the  stimulated 
Raman  effect  are  analyzed,  and  some  others  arising  from  transverse 
Inhomogeneity  are  identified. 
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INTRODUCTION 


The  recent  successful  operation  of  the  "free  electron  laser"  by 
the  group  at  Stanford  and  the  availability  of  high  power  electron  beans  has 
stimulated  a  great  deal  of  interest  in  the  use  of  the  "free  electron  laser" 
(FEL)  to  produce  a  high  power  tunable  laser  beam.l  When  an  electron  tra¬ 
vels  through  a  periodic  transverse  magnetic  field,  the  electron  is  given  a 
transverse  velocity  which  allows  it  to  either  receive  or  give  energy  to  the 
transverse  electric  field  of  a  plane  electromagnetic  wave.  If  the  longitu¬ 
dinal  velocity  of  the  electron  is  such  that  the  electron  slips  behind  the 

radiation  wave  by  one  radiation  wave  length  while  traveling  a  distance  of 

one  magnetic  field  period,  the  transverse  velocity  of  the  electron  remains 
in  phase  with  the  electric  field  of  the  radiation.  Such  electrons  will 
continue  to  have  their  energy  Increased  or  decreased,  depending  on  their 
phase  relative  to  the  radiation  field.  Phillips  investigated  such  a 
device,  called  the  Ubitron,  and  used  it  to  produce  microwaves. ^  The  FEL 
experiment  by  the  Stanford  group^  used  a  helical  wiggler  magnet  with  a 
constant  magnetic  field  and  periodicity  and  produced  optical  radiation  in  a 

similar  manner  as  the  Ubitron.  Most  of  the  analysis  of  the  FEL,  as  well  as 

the  methods  envisioned  for  FEL  operation,  Btarted  with  the  assumption  that 
the  magnetic  field  and  period  of  the  wiggler  are  constant. 5  a  brief 
review  of  the  properties  and  limitations  of  such  a  device  is  presented  in 
Section  3  of  this  paper. 
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The  maximum  fractional  energy  that  can  be  extracted  from  a  con¬ 
stant  parameter  FEL  is  equal  to  1/2N  where  N  is  the  number  of  wiggler 
periods  while  the  resulting  energy  spread  beam  of  the  electrons  emerging 
from  the  FEL  is  comparable  to  or  greater  than  the  mean  .energy  loss.  This 
places  severe  limitations  on  the  overall  gain,  efficiency,  and  modes  of  c 
operation  for  the  constant  parameter  FEL.  While  there  are  many  schemes  to 
overcome  some  of  the  limitations  of  the  constant  parameter  FEL,  such  as  the 
use  of  transverse  gradients  in  the  magnetic  field  of  the  wiggler®  and  the 
coupling  of  the  FEL  with  an  isochronous  storage  ring,^  the  use  of  a 
variable  parameter  wiggler  may  be  an  even  more  attractive  method  to  remove 
these  limitations. 

It  is  the  purpose  of  this  paper  to  discuss  the  operation  of  a 
variable  parameter  FEL  and  to  illustrate  many  of  the  possible  operational 
modes  that  may  be  used  to  overcome  these  limitations  of  the  constant  para¬ 
meter  FEL.  For  the  case  where  the  wiggler  parameters  such  as  magnetic 
field  and  wavelength  are  variable,  the  possible  operating  modes  envisioned 
for  an  FEL  are  considerably  different  than  those  for  a  constant  parameter 
wiggler;  for  example,  it  becomes  possible  to  extract  a  reasonable  fraction 
of  the  electron's  energy  with  an  emerging  energy  spread  a  great  deal  less 
than  the  mean  energy  loss.  The  method  of  describing  the  electrons'  motion 
in  this  paper  is  similar  to  the  treatment  used  by  the  accelerator  physicist 
in  describing  the  motion  of  charged  particles  in  a  radio  frequency 
accelerating  system.®*^  We  assume  from  the  beginning  that  the  magnetic 
field  and  period  vary  along  the  wiggler  and  derive  the  equations  of  motion 
for  an  electron  passing  through  the  FEL,  from  which  we  can  obtain  the 
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energy  lost  or  gained  by  the  electron.  The  change  In  the  radiation  field 
is  taken  to  be  given  by  the  sum  of  the  energies  lost  by  all  of  the 
electrons  which  pass  through  the  FEL. 

In  Section  2.1  we  present  the  equations  that  define  the  energy  and 
relative  phase  of  a  resonant  or  synchronous  electron  in  terms  of  the  wig- 
gler  magnetic  field,  the  wiggler  period,  the  optical  field,  and  the  optical 
wave  numbers.  In  Section  2.2  we  discuss  the  motion  of  the  electrons  about 
the  synchronous  energy,  using  the  Hamiltonian  formulation  similar  to  that 
used  in  the  treatment  of  charged  particle  acceleration  with  radio  frequency 
fields.  The  summary  of  the  "bucket"  or  maximum  stable  phase  trajectories 
for  the  electrons  is  presented  in  Section  2.3.  The  results  of  Section  2 
are  used  in  Section  3.0  to  briefly  discuss  the  well-known  properties  of  the 
constant  parameter  wiggler,  both  as  an  illustration  of  the  use  of  the 
Hamiltonian  formulation  and  to  obtain  the  results  in  terms  of  notation  of 
this  paper. 

In  Section  4  the  use  of  a  variable  parameter  wiggler  to  trap  a 
significant  fraction  of  the  electron  beam  in  a  decelerating  bucket  is 
discussed.*®  The  wiggler  parameters  are  varied  so  as  to  gradually  reduce 
the  resonant  energy  of  the  trapped  electrons,  which  results  in  a  value  for 
the  transfer  of  energy  from  the  trapped  electrons  to  the  radiation  approxi¬ 
mately  equal  to  the  reduction  in  the  resonant  energy.  This  reduction  in 
the  resonant  energy  can  be  much  larger  than  the  energy  that  can  be  trans¬ 
ferred  in  the  constant  parameter  wiggler.  The  operation  of  a  laser 
oscillator  which  utilizes  this  scheme  is  discussed  in  Section  4.1  where  the 


amplification  of  the  radiation  field  in  passing  through  the  wiggler  just 
compensates  for  the  reflection  loss.  Amplifier  designs  for  this  trapping 
and  deceleration  mode  of  FEL  operation  are  discussed  in  Section  4.2.  While 
the  trapped  electrons  emerge  with  an  energy  spread  comparable  to  that  of 
the  constant  parameter  FEL,  the  mean  energy  loss  of  electrons  is  an  order 
of  magnitude  larger  than  that  from  the  constant  parameter  FEL. 

In  Section  5  a  method  of  adiabatically  capturing,  deceleration, 
and  decapture  of  the  electrons  is  discussed.  The  purpose  of  this  scheme  is 
to  achieve  a  small  energy  spread  for  the  emerging  beam  while  at  the  same 
time  to  obtain  a  significant  reduction  of  the  average  electron  energy. 

Such  a  scheme  would  be  useful  if  the  same  electron  beam  was  reaccelerated 
with  an  external  source  and  then  passed  through  the  FEL  again. 

The  last  method  of  extracting  energy  from  the  electrons,  discussed 
in  Section  6,  is  called  phase  area  displacement,  and  is  one  in  which  the 
initial  resonant  energy  is  below  the  energy  of  the  injected  electrons. ^ 

The  wiggler  parameters  are  varied  such  that  the  resonant  energy  is  in¬ 
creased  to  a  final  value  above  the  energy  of  the  electrons,  a  process  which 
can  result  in  a  decrease  in  the  mean  energy  of  the  electrons  while 
producing  a  rather  modest  increase  in  the  energy  spread  of  the  beam.  While 
the  average  energy  lost  by  this  process  for  mono-energetic  electrons  is  of 
the  same  order  as  the  energy  lost  from  a  constant  parameter  FEL,  it  has  the 
advantage  of  being  less  sensitive  to  the  initial  energy  spread  of  the 
electrons  and  of  producing  a  lower  energy  spread  in  the  emerging  beam. 
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There  is  a  possibility  of  an  Instability  that  can  arise  from  the 
interaction  between  a  radiation  signal  at  a  frequency  shifted  from  the  main 
radiation  frequency  by  an  amount  equal  to  the  oscillation  frequency  of  the 
electrons  in  the  pondermotive  potential  well;  this  instability  is 
investigated  in  Section  7.  The  results  show  that  the  instability  that  must 
build  up  from  noise  is  probably  not  important  for  the  amplifier,  but  in  the 
design  of  an  oscillator  some  type  of  damping  will  probably  need  to  be 
included. 


Some  related  considerations  are  presented  in  the  Appendices. 
Appendix  A  gives  a  more  general  derivation  of  the  Madey  Theorem1^  relating 
small  signal  energy  gain  and  spread  with  application  to  the  variable 
wiggler.  Appendix  B  presents  some  mathematical  details  relating  to  Chapter 
VI.  Appexdix  C  is  concerned  with  the  requirements  on  e-beam  emittance. 
Appendix  D  discusses  the  justification  at  high  Y  for  neglecting 
electrostatic  effects.  Appendix  E  discusses  a  different  variable  wiggler 
concept,  the  Klystron, ^  which  is  related  to  ff-pulse  designs,1^1  in  which 
one  wiggler  is  used  for  prebunching  and  a  second  for  energy  extraction. 
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PRECEDING  FACS  BLANK. NOT  FILMED 


2.0  EQUATIONS  OF  MOTION 

The  equations  of  motion  of  an  individual  electron  are  determined 

by  a  Hamiltonian,  tfo<*.  Px.  F.  Py  ?z*  which  8atisfies  a  nodified 
Hamiltonian  variational  principle 

c2 

6  J  (px*  +  P^y  +  P^2  -  HQJdt  ■  0  (2.1) 

*1 

Throughout  this  derivation  we  shall  assume  that  the  electrostatic 
interaction  between  the  electrons  can  be  neglected,  deferring  until  later 
the  determination  of  the  circumstances  in  which  this  assumption  is 
justified.  With  this  assumption  we  can  write 

Hq( r,  P,  t)  -  c  J  m2c2  +  [p  -  (e/c)  A(r ,  t) J2  (2.2) 

where  the  total  energy,  E  ,  of  the  electron  Is  to  be  identified  with  the 
Hamiltonian,  i.e.  E  -  H0  •  We  expect  many  of  the  various  parameters  to 
vary  with  the  longitudinal  coordinate  z  ;  for  this  reason  it  is  useful  to 
identify  z  as  the  independent  variable  rather  than  the  time  t  •  We 
rewrite  Eq.  (2.1)  as 

z2 

6  [  [p  x'  +  P  y'  +  (-E)t'  -  (-P  ) ]dz  -  0  (2.3) 

J  x  y  z 

Z1 

where  we  define  the  derivatives  with  respect  to  z  with  a  prime,  i.e.  x' 
■  dx/dz  ,  etc. 
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as  a  new 


It  is  convenient  to  designate  the  quantity  (-Pz) 
Hamiltonian  H ^  such  that 


Hl(x’  Px*  y>  Py*  (_E>*  tf  z)  “  ”P2 


(2.4) 


where  (-E)  plays  the  role  of  the  momentum  conjugate  to  t  and  z  is  the 
independent  variable.  The  solution  of  Eq.  (2.2)  for  yields 


(E2/c2  -  m2c2]  -  (p  -  e/c  A  )2  -  (p  -  e/c  A Y 
x  x  y  y 


1/2 


-  e/c  A 


(2.5) 


For  the  vector  potential  we  shall  write 


A  -  A  +  A 


(2.6) 


where  A  refers  to  the  vector  potential  of  the  wlggler  field  and  4  to 
"  s 

that  of  the  signal  (or  optical)  field.  We  assume  these  potentials  to  have 
the  special  forms 


A  -  -  A  (z) 
-w  w 


X  cos 


z  z 

t/kw(zl)dzlJ  +  y  sin(  Jv.^(z^)dz  ^  ) 


(2.7) 


As(z) 


X  COS 


</ 


ka(zi)dzi 


-  ■ 


y  sin( / kfi(z1)dz j 
o 


(2.8) 
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Equation  (2.7)  corresponds  to  a  circularly  polarized  static  magnetic  field 

of  the  sort  used  in  the  original  Stanford  FEL  experiments.  The  assumption 

of  no  transverse  variation  can  only  be  an  approximation  in  practice,  but  we 

assume  that  the  electrons  are  confined  to  a  region  in  which  transverse 

variation  can  be  neglected.  The  use  of  a  non  static  or  linearly  polarized 

wlggler  field  makes  no  essential  difference  In  the  theory,  but  the 

assumption  of  the  form  (2.7)  allows  some  simplification  of  the 

exposition.  A ^  and  ky  are  both  taken  to  be  z  dependent  to  allow  for 

the  use  of  specialized  magnet  designs.  Ue  assume  A  ,  k  ,  A  ,  k  to  be 

w  w  s  s 

slowly  varying  in  z  in  the  sense  that  it  will  generally  be  permissible  to 

neglect  A'  -  dA  /dz  as  compared  to  k  A  ,  etc.  The  signal  field  has 
w  w  r  w  w 

also  been  chosen  to  be  circularly  polarized  (in  the  direction  that  is 
driven  by  the  FEL  amplification  process). 

The  assumption  that  A  is  independent  of  the  transverse 

A  A 

coordinate  implies  that  p,  «  Px+Py  is  a  constant  of  the  motion.  We 

'“i.  x  y 

shall,  throughout  most  of  this  paper,  assume  this  constant  to  be  zero,  an 
assumption  which  corresponds  to  the  neglect  of  transverse  velocities  in  the 
incoming  beam  before  it  enters  the  wiggler.  Some  aspects  of  the  effect  of 
non  zero  P^  on  operation  are  discussed  in  Appendix  C  and  applied  to  a 
discussion  of  the  relation  of  the  characteristics  of  realizable  electron 
beams  to  other  design  parameters. 

Taking  account  of  the  above  assumptions  we  find  that  the  equations 
of  motion  for  the  energy  and  time  can  be  obtained  from  the  following 
Hamiltonian 

9 


w 


((-E) ,  t,  z )  m  -mC 


-U2(z) 


T2  *  2V, 


C08  <P(Z,  t) 


jl/2 


(2.9) 


where  we  have  Introduced  the  dimensionless  vector  potentials 

2  2  6 
a  -  e  A/mc  ,  (me  /e  -  1706  gauss-cm  ■  .511  *  10  Volts) 


P2(z)  »  1  +  a2  (z)  +  a2(z)  , 

z 

^  ■  f  (k  +  k  )dz.  -  u  t  and 

J  W  8  1  3 

O 

Y  -  E/mc2 


(2.10) 


(2.11) 


(2.12) 


With  the  special  assumptions  that  we  have  made,  the  problem  has  been 
reduced  to  a  straightforward  one  of  one  dimensional  relativistic  motion 
with  specified  external  forces.  The  quantity,  mV  may  be  thought  of  as  an 
effective  electron  mass,  increased  above  the  rest  mass  by  the  transverse 
kinetic  energy. 


From  Hamilton's  equations  (in  the  context  of  Eq.  2.9), 

3H^/3(-E)  -  t'  ,  3^/3t  -  -(-£')  ,  and  the  definitions  of  ip  and  y  Eqs. 
(2.11)  and  (2.12)  we  obtain 


r 


(k  +  k  )  + 

W  8 


U), 


me 


■s 

2  3y 


(2.13) 


and 


Y' 


Wg  3  Hj 

- ITT 

me 


(2.14) 
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Strong  beam-wave  interaction  occurs  when  the  beam  velocity  is  approximately 
matched  to  that  of  the  bunching  wave,  in  which  case  ♦  is  approximately 
constant.  While  we  expect  to  allow  the  various  parameters  to  vary  with 
z  it  will  still  be  the  case  that  the  orbits  of  principal  Interest  will  be 
those  which  remain  in  approximate  velocity  resonance  and  for  which  ♦  is 
slowly  varying.  For  this  reason  it  proves  to  be  more  convenient  to  use 
4*  and  Y  as  (canonical)  coordinates  rather  than  t  and  (-E)  with  z  still 
the  independent  coordinate.  This  leads  us  to  use  the  following  Hamiltonian 
to  describe  the  motion  in  the  Y  coordinates 


H 

2 


Y(k  +  k  ) 
w  s 


+  2a  a  cos 
w  s 


* 


1/2 


(2.15) 


The  equations  of  motion  for  ♦  and  Y  which  follow  from  this  new 
Hamiltonian  are 


V 


(k  +  k  )  - 

W  8 


“s 


J Y^  -  +  2a  a  cos 

V  w  s 


(2.16) 


Y  -  - 


u 

T-  Vs  8ln  * 

. /  Y^  -  +  2a  a  cos  ♦ 

*v  w  s 


(2.17) 


which  of  course  are  the  same  equations  that  are  obtained  by  the  use  of  Eqs. 
(2.13)  and  (2.14)  with  the  old  Hamiltonian  H^(-E,  t,  z)  • 


2  2 

For  the  applications  to  be  discussed  in  this  paper  y  >>  W  , 

Y^  »  *w*s  and  *w  >>  *g  •  With  these  approximations  the  Hamiltonian  may 
be  approximated  by 
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(2.18) 


H3  "  Y(kw  +  6ks)  +  773  <W 


2a  a  cos  v) 

V  8 


with 


M2  *  1  +  a2  and  <5k  -  k  -  ^0.  . 

w  s  s  c 


(2.19) 


Evidently  the  term 


U  - 

P 


a  a  u) 
w  s  s 

Yc 


k  )dz, - 
w  1 


V 


(2.20) 


plays  the  role  of  a  potential  (the  "ponderomotive  potential").  It  may  be 
thought  of  as  providing  a  bunching  force  which  moves  at  the  resonant 
velocity. 


v 

r 


k  +  k 


(2.21) 


The  equations  of  motion  that  correspond  to  the  approximate  Hamiltonian 
H3  are 


<»'  "  (k  +  <5k  )  - — 1~  (w2  -  2a  a  cos  i|>)  , 

w  s  w  s 


(2.22) 


and 


Y'  -  - 


usasaw 

YC 


sin  v 


(2.23) 
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The  corresponding  equations  for  the  particle  energy  and  phase  nay  be 
written  as 


u  11  -  2a  a  cos 
s  w  s 

2c  7k  +~Sk  -  r) 


co 

s 


2c(k  + 
w 


s 


*  > 


(2.24) 


and 

2(k  +  6k  -  *')a 

-  (k'  +  6k') - - - 1 - —[a'  +  (2k  +  2«k  -  *')a  sin  ♦] 

w  s  ^2  1  w  w  s  s 

(2.25) 

where  we  have  used  the  approximation  that  a  »  a 

V  8 

The  z  dependence  of  a^  and  ky  are  determined  by  the  physical 
design  of  the  FEL  hardware,  and  in  subsequent  sections  we  will  discuss  the 
design  considerations  upon  which  these  specifications  may  be  based.  The 
optical  field  amplitude,  and  to  some  extent  the  phase,  depend  explicitly 
upon  the  FEL  amplification  process  and  will  be  affected  by  such  additional 
properties  as  the  design  of  the  optical  cavity  and  the  design  and  magnitude 
of  the  electron  beam.  In  keeping  with  our  assumption  of  the  absence  of 
variation  in  the  transverse  direction  we  shall  assume  that  the  electron 
beam  and  optical  beam  overlap  perfectly  and  that  the  electron  beam  is  also 
uniform  in  electron  density.  Under  these  conditions  the  variation  of  ag 
is  determined  by  energy  conservation.  Thus  we  write 
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x  B 


-s  z 


*  s,u 

8  0 


-  c  k2 
T*  s 


2  2 


2 

■  n  me 
e 


(<Y(0)>  -  <Y  (z)  >  J 


or 

2  2  WD 

as(2)  ’  aa(0)  +-7  -  <Y (z ) > j  (2.26) 

s 

where  <  >  implies  averages  over  initial  distributions  of  energy  and  phase 
and  <j2  .  ATTnee2/m  is  the  plasma  frequency.  n  is  the  electron 
density  measured  in  the  laboratory  frame,  and  we  have  assumed  v  »  c  • 


An  alternative  derivation  of  Eq.  (2.26)  along  with  a  determination 
of  6ka  can  be  obtained  from  the  Maxwell  equation 


f*. .  i  3X 
<7  7  3t^ 


(2.27) 


Substitution  of  the  form  (2.8)  yields 


T~3  6ksAsel  *  2  ^  A;e2  "  £  Jj. 


(2.28) 


where  and  e^  are  defined  by 


z  x 

!1  “  [x  coa  </  k8<z1)  dz!  -  «,t>  -  y  sin  {f  k8(Zj)  dZj  -  u>gt) 


*i  -  kse2 


(2.29) 
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and  we  have  neglected  derivatives  of  A  '  and  gi.  compared  to 

S  *8 

a  a  a 

Because  e^  is  a  unit  vector  and  e^*  vanishes  we  find 


w8/c  • 


2iufi 


<Sk  A 
s  s 


(2.30) 


In  order  that  the  form  originally  chosen  for  A  be  strictly  correct  it  is 

~s 

necessary  that  the  right  hand  side  of  Eq.  (2.30)  be  time  independent.  In 

actual  fact  the  generated  by  the  assumed  form  for  As  will  be  a 

periodic  function  of  time  with  period  2*/w  .  We  shall  eliminate  the 

s 

oscillating  terms  by  carrying  out  a  long  time  average  of  •  e^  •  Such 
a  procedure  is  valid  ^  provided  6kg  «  kw 

For  a  single  electron  we  have 


*  ej  -  ev±  •  e^6(x  -  xQ(t)  )<5(y  -  yQ(t)  )<5(z  -  z(t) )  (2.31) 


where 


v,  - 


V 

~P  1 


-eA  -eAw 

m  . 

Ymc  Yrac 


(2.32) 


for  A  <<  A  .  Hence, 
s  w 

e?A  (z) 

(il  *  *1  Single  particle  THc  008  *  "  xo(t))6[y  “  yo(t))6(z  "  zo(t) ] 

(2.33) 
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We  time  average  by  integrating  over  time  from  -T/2  to  T/2  to  obtain 
-  e^A  (z) 

’  *1  'single  particle  '  C0*  '>(2)  5(x  '  V2>)  ‘(y  ‘  V2,) 

Z 


for  -T/2  <  t(z)  <  T/2 


0  otherwise. 


(2 . 34) 


Summing  over  all  electrons,  and  averaging  over  the  beam  cross 
section  to  eliminate  the  transverse  6  functions  we  obtain 


(i,  •  e,) 


2 

n  ezA 
e  w 


1  1  time  average  me 


<cos  ii>/Y> 


(2.35) 


where  again  <  >  means  average  over  initial  ^  and  energy.  Substitution 
of  Eq.  (2.35)  with  Eq.  (2.30)  yields 


6k 


w  ^  a 
p _ w 


s  2cu  a 
s  s 


<cos  <|>/Y> 


(2.36) 


Again  from  Eq.  (2.28)  we  find 


2  _s  A.  .  .  41  ; 

c  s  c  2 


and  proceeding  in  a  similar  manner  we  find 


u> 

a'  »  —  a  <sin  <l>/ Y> 

s  2w  c  w 
8 


(2.37) 
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To  compare  Eq.  (2*37)  with  Eq.  (2 .26)  we  note  that 


dY 

Sz 


— z  ■ar’ 

me 


a  a  (i) 

— 8CW--  (sin  i|>/Y) 


so  that 


d<Y> 
d  z 


asawus 


<sin  'J'/Y> 


(2.38) 


The  equivalence  of  Eq.  (2.37)  and  Eq.  (2.38)  to  Eq.  (2.26)  is  apparent. 
From  Eq.  (2.36)  and  Eq.  (2.37)  we  have 


6k  <cos  *IVy> 

s  “  a  <8in  <i*/Y> 
8 


(2.39) 


so  that  the  assumption  a' /a  «  k  made  at  the  outset  would  appear  to 

8  8  W 

typically  imply  6kg  «  kw  as  well.  Since  we  expect  the  particles  to 
bunch  in  the  0  <  ^  <  x/2  range,  <cos  ♦/Y>  will  be  positive.  This 
implies  a  tendency  of  the  electron  beam  to  trap  the  optical  beam  and  hence 
to  counter  to  some  extent  the  effects  of  diffraction. 


The  equations  of  motion  derived  above  (Eqs.  (2.24)  and  (2.25)  )  are 
valid  (under  the  special  assumptions  that  were  made)  for  every  electron, 
and  given  the  initial  conditions  <Ko)  and  Y(o)  for  every  electron 
along  with  the  expressions  for  ky,  Sy,  a8,  and  6  kg  ,  these  equations 
may  be  integrated  to  yield  the  values  of  ♦  and  Y  as  functions  of  the 
longitudinal  position  z  .  Of  course,  for  a  high  gain  FEL  it  is  necessary 
to  make  sure  that  the  signal  field  a8  used  in  the  equations  of  motion  is 


also  a  self  consistent  solution  to  Eqs.  (2.26)  and  (2.39);  by  using  a 
computer  one  can  solve  the  equations  of  motion  for  a  large  number  of 
particles  and  obtain  the  final  coordinates  In  phase  space  for  every 
particle  as  well  as  the  final  signal  field  and  phase.  In  principle  this 
procedure  could  be  used  to  choose  the  optimum  functional  form  for  k^z) 
and  aw(z)  .  It  is  possible  to  gain  physical  insight  Into  the  content  of 
these  equations  without  the  need  of  Integrating  the  equations  of  motion  for 
a  large  number  of  initial  conditions  by  referring  the  electrons  energy 
variable  Y  to  the  synchronous  or  resonant  value  Yf  •  This  approach  is 

extremely  useful  in  determining  the  functions  kw(z)  and  aw(z)  to  be 
used  in  different  modes  of  operation  for  the  FEL. 


2.1  Definition  of  Synchronous  Energy  and  Phase 

The  synchronous  energy,  Yr  ,  and  phase,  ^  ,  are  defined  by 


„  U)  u* 

\  "  2"(k  +  5k  )c 

W  8 


(2.40) 


and 


d\ 

ST 


a  a  w 
w  s  s 

Y  c 

r 


sin  ♦_ 


(2.41) 


If  Eq.  (2.41)  yields  a  z  independent  value  for  ,  then  substitution  in 
Eqs.  (2.24  &  2.25)  shows  that  a  particle  with  values  of  Y  and  V  equal 
to  Y^  and  4*r  respectively  will  retain  these  values  throughout  its 
motion  through  the  wlggler.  Such  a  particle  may  be  thought  of  as  a 
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synchronous  particle.  More  generally,  as  will  become  apparent  below,  if 
Y  -  y  ,  -  i>T  are  initially  small,  they  tend  to  remain  small. 

It  is  possible  to  look  at  Eqs.  (2.40)  and  (2.41)  as  definitions 

of  Yr  and  *r  assuming  that  kw,  aw,  as  and  6kg  are  known  functions 

of  z  .  However,  it  is  also  possible  to  consider  these  as  design  equations 

where  the  wiggler  functions  and  are  to  be  determined  to  achieve 

the  desired  functions  Y  ,  4*  and  a_  • 

r  r  3 

Much  of  the  remainder  of  the  paper  will  deal  with  different  ways 
to  choose  and  Yr  and  hence  the  wiggler  parameters  for  various 

operating  modes  of  the  FEL.  We  first  must  study  the  motion  of  electrons 
with  phase  coordinates  different  from  the  synchronous  values  and 
demonstrate  that  for  sufficiently  small  deviations  the  electrons  will 
perform  stable  oscillations  about  the  synchronous  values.  It  must  be  noted 
that  one  is  restricted  in  the  choice  for  the  functions  Yf  and  ^  ;  for 
example  ♦  is  not  defined  if 


dY 

w  a  a 

r 

> 

s  w  s 

dz 

CY 

r 


2.2  Motion  About  the  Synchronous  Energy 

In  this  section  we  study  the  motion  of  the  electrons  about  the 
synchronous  energy  by  writing 

Y  -  $Y  +  Y  •  (2.42) 

r 
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We  regard  6y  as  the  new  canonical  momentum  with  a  new  Hamiltonian  ^ 
given  by 


J 


<^r 


<5y,  z)  -  H3(4S  Yr  +  «Y,  z)  +  *  dz 


(2.43) 


With  the  definitions  for  Y  and  ^  as  chosen  in  Eqs.  (2.40)and  (2.41) 

r  r 


and  with  the  assumption  of  Sy  «  Yr  the  new  Hamiltonian  ^  f  which  is 
derived  from  the  approximate  Hamiltonian  H 3  Eq.  (2.18),  is 


K  +  6ks 


2  u>  a^ 

($Y) - - —  (cos  '{»+'{'  sin  ^  ) 

c  Y  r 


(2.44) 


This  form  for  the  Hamiltonian  and  the  resulting  equations  of  motion  are 
very  familiar  to  accelerator  physicists.  The  following  analysis  is  quite 
similar  to  their  treatment  of  RF  acceleration  in  standing  wave  linacs.  In 
writing  Eq.  (2.44)  we  have  omitted  terms  which  are  functions  only  of  the 
independent  variable  z  ,  since  they  have  no  effect  upon  the  equations  of 


motion,  and  it  has  been  assumed  that  a  value  of  exists,  i.e.  that 


He 

„  a)sawas 

dz 

Yc 

(2.45) 


H4  has  the  form  of  a  nonrelativistic  single  particle  Hamiltonian  with  z 
dependent  "mass"  and  a  z  dependent  potential  function,  F,  given  by 


F(\|»)  *»  -C(cos  +  'J'  sin  4^) 


(2.46) 


A 
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with 


C(z) 


U)sasaw 

V 


(2.47) 


We  assume  C  >  0  and  -ir/2  <  4*  <  ti/2  ,  so  that  corresponds  to  a 

r  r 

minimum  of  the  potential  and  (it  sgn  ♦  -  ij^)  to  a  maximum.  There  are,  of 
course,  a  succession  of  minima  at  ♦  ■  ^  +  2*  n  and  maxima  at 
♦  *  it  -  +  2ti  n  ,  and  one  may  conveniently  think  of  the  successive 

troughs  as  buckets  in  which  particles  may  become  trapped.  F^)  is 
illustrated  in  Fig.  2.1. 


Of  particular  interest  are  design  parameters  chosen  so  that  the 
behavior  changes  adiabatically  with  z  .  In  that  case  electrons  trace  out 
trajectories  in  the  Y,  5Y  phase  plane  given  by 


6Y(H,  *,  z) 


+ 


[H  -  F('j')  ] 


(2.48) 


and  illustrated  for  0  in  Fig.  2.2  at  z  ■>  0  .  The  case 

0  corresponds  to  dYf/dz  <  0  and  therefore  represents  decelerating 
buckets.  ”  (0)  is  determined  by  the  initial  values  of  ^  and  •  The 

closed  orbits  correspond  to  the  particles  trapped  in  the  buckets.  If  the 
change  of  parameters  with  z  is  adiabatic,  then  H  (z)  is  determined  by 
the  requirement  that  the  area  of  the  closed  phase  curve,  given  by 


J  - 


/ 


6ird>|* 


(2.49) 
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FIGURE  2.2 


TRAJECTORIES  IN  THE  i,  6>  PHASE  PLANE  FOR  f 


remain  constant  as  z  changes*  For  the  unbound  orbits,  the  representation 
Is  valid  only  over  a  range  of  z  for  which  the  change  In  parameters  Is 
small.  For  small  oscillations  about  4*r  one  can  expand  F(^)  about  ^  The 
motion  for  these  orbits  Is  harmonic  with  period  of  oscillation 


Z 


HU 


(k  +  5k  ) 
w  s 


J 


a  a  cos 
s  w 


pi 

_ w _ 

2  /  a  a  cos 

s  w  r 


2n/kw)  . 


(2.50) 


The  period  for  weakly  trapped  particles,  with  large  excursion  In  ^  ,  Is  of 

course  larger.  The  parameter  variation  may  be  considered  to  be  adiabatic 

If  It  is  small  over  a  distance  of  the  order  of  the  period.16  In  order  that 

a  particle  be  trapped  In  a  bucket  it  Is  necessary  that  |  <5 Y |  <  <$Y  with 

max 


oY  o  I  “ 

ym-a~X'  *  -y  Javaa  lcos  “  (7  s8n  ~  *r)  sin  <1^]  (2.51) 

If  one  wishes  to  avoid  the  approximations  associated  with  the  use  of 
Eq.  (2.44)  It  is  straightforward  to  base  the  definition  of  Y^  ,  and 
J  upon  Eqs.  (2.15  &  2.43). 


2.3  Summary  of  Bucket  Parameters 

Because  the  pattern  of  the  stable  trajectories  traced  out  by  the 
electrons  in  the  t,  6y  phase  plane  repeat  in  ♦  at  Intervals  of  2*  we  only 
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.1 


will  discuss  the  phase  notion  in  the  interval  for  -  r  <  <  it  .  The 

maximum  stable  phase  curve  or  bucket  is  shown  in  Fig.  2.3  and  corresponds 
to  the  trajectory  given  by  Eq.(2.48)  for 


H  -  H  -  C  [cos  ip  -  (nsgn  ^  -  ip  )  sin  ip  ]  (2.52) 

m  r  r  r  r 


wsasaw 


(2.53) 


The  bucket  intercepts  of  the  ip  axis,  designated  by  ip^  and  ip2  in  Fig.  2.3, 
are  given  by 


ip2  =  "sgn  ipr  -  ip^ 


(2.54) 


cos  ip,  +  ip,  sin  ip  ■  cos  ip„  +  ip„  sin  ^ 

1  1  r  2  2  r 


(2.55) 


The  values  of  ip^  and  are  shown  as  a  function  of  in  Fig.  2.4.  The 

maximum  bucket  height  given  by  Eq.  (2.51)  can  be  written  as 


r  (.  > 


(2.56) 


r(  +r )  ■  Jcos  ipr  -  (^  sgn  ipr  -  ipr)  sin  ipr 


(2.57) 


T  (ip^)  is  shown  in  Fig.  2.5. 
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3661 A2 


FIGURE  2.3 


STABLE  PHASE  PLANE  TRAJECTORIES 


The  area  of  the  bucket  shown  in  Fig.  2.3  nay  be  obtained  from  Eq. 


(2.49)  and  is  given  by 


, _  f2  r  ■  12 

2  p  V  2asaw  J  Jcos  +  cos  ^  -  [n  -  Vr  -  Vj  sin  vj  d  v 


(2.58) 


16  Yr  / -  .  . 

_ _ _ J  a  a  a(v  ) 

p-  v  s  w  r 


(2.59) 


where 


(vr) 


ms/2 

TT 


P  r 

#  cos 

L 


V  + 
r 


cos  \|/  -  (n  -  ^  -  ^)  sin  ij» 


(2.60) 


is  the  moving  bucket  area  and  is  plotted  in  Fig.  2.6.  Note  that  a  equals 
one  for  «  0  and  decreases  to  zero  for  $tm  +it/2  .  ^Eqs.  (2.58)  and 

(2*60)  are  for  Vr  >  0  ,  For  <  0  ,  use  a(^)  -  a(-  Vy)  . ] 
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3.0  CONSTANT  PARAMETER  WIGGLER 

O 

The  original  mode  of  FEL  operation  demonstrated-*  by  Madey  and  co¬ 
workers  utilized  a  wiggler  with  a  fixed  wave  number  and  field 

amplitude^  aw  .  For  such  a  case  the  operational  mode  is  one  in  which 


* 


r 
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where  the  fact  that  6k  «  k  has  been  used. 
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The  Hamiltonian^  which  describes  the  motion  of  the  electrons  is 


given  by 
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For  this  mode  of  operation  the  buckets  are  stationary  or  non-accelerating 
with  a  half-height 


and  the  length  of  one  synchrotron  period 
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The  motion  of  the  electrons  in  such  a  bucket  is  one  in  which  some  of  the 
electrons  gain  energy  while  others  lose  energy  in  a  manner  that  depends 
upon  the  total  length  of  the  FEL  as  well  as  the  initial  phase  and  energy 
of  the  electron.  If  the  average  injection  energy  equals  the  resonance 
energy,  the  average  change  in  the  electron  energy  will  be  zero,  provided 
that  the  initial  phase  distribution  is  uniform.  Whi  i.e  there  are  modes  of 
operation  for  the  constant  parameter  wiggler  which  assume  a  non-uniform 
phase  distribution  for  the  injected  electrons,  we  will  consider  only  the 
case  where  the  initial  phase  distribution  is  unifcrm. 


Electrons  that  have  an  initial  energy  different  from  the  resonant 
energy  by  an  amount  large  compared  to  the  maxiumum  bucket  height  will  not 
be  appreciably  affected  by  the  FEL  and  hence  will  have  only  a  small  change 
in  their  energy  (they  will  not  be  resonant  with  the  optical  wave),  while 
electrons  that  are  too  near  the  resonant  energy  also  will,  on  the  average, 
have  only  a  small  energy  change.  The  electrons  with  an  initial  energy 


difference,  from  resonance,  (6y)^  will  have  the  largest  energy  change  in 
traveling  through  the  FEL  when 


(6Y)  ~  ( fiy  ) 

1  r  m 


Clearly,  if  a  significant  fraction  of  the  electrons  are  to  lose 
energy,  it  is  necessary  for  the  initial  energy  spread  in  the  bean  to  be 
less  than  the  maximum  bucket  height.  Consider  a  mono-energetic  beam  of 
electrons  with  an  initial  energy  difference  (<5y)^  >0  as  shown  in  Fig. 
3.1a.  After  a  distance  L  which  corresponds  to  approximately  one-half  of 
a  synchrotron  oscillation  many  particles  will  have  their  energy  shifted 
below  the  resonant  energy  as  shown  in  Fig.  3.1b.  For  longer  distances 
these  electrons  will  continue  to  oscillate  about  the  resonant  energy  and 
the  average  energy  will  begin  to  increase  back  to  the  initial  energy  as 
shown  in  Fig.  3.1c. 


The  performance  of  a  FEL  operating  with  an  electron  beam  of  an 
initial  uniform  phase  distribution  and  a  constant  parameter  wiggler  will 
have  the  following  properties.  First,  for  reasonable  performance,  the 
energy  spread  of  the  beam  must  not  be  significantly  larger  than  the  bucket 
height,  secondly,  the  maximum  energy  loss  by  the  beam  will  be  of  the  order 
of  the  bucket  height  and  thirdly  the  optimum  length  of  the  FEL  is  of  the 
order  of  one-half  of  a  synchrotron  period.  From  Eqs.  (3.5  &  3.6)  we  see 
that  the  maximum  allowable  energy  spread,  maximum  energy  loss,  and  optimum 
FEL  length  are  related  by 
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FIGURE  3.1 


PHASE  SPACE  PLOTS  FOR  MONO-ENERGETIC  ELECTRONS  IN  A 
CONSTANT  PARAMETER  WIGGLER. 

a)  at  z  -  0 

b)  L  corresponds  to  ~  1/2  synchrotron  oscillation 

c)  L  is  5/3  the  value  shown  in  b. 
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where  N  is  the  number  of  wiggler  periods. 


Another  important  characteristic  of  the  constant  parameter  wiggler 
is  the  fact  that,  due  to  the  dependence  of  final  energy  upon  initial  phase, 
electrons  which  enter  at  a  particular  energy  emerge  with  an  energy  spread 
that  is  of  the  same  order  (typically  larger)  as  the  mean  energy  loss.  This 
property  is  qualitatively  apparent  in  Fig.  3.1b.  As  a  result  the  energy 
spread  in  the  emerging  electrons  tends  to  be  increased,  over  that  of  the 
entering  electrons,  by  an  amount  comparable  to  the  mean  energy  loss.  The 
above  refers  to  the  situation  in  which  the  signal  amplitude  is  such  that 
the  synchrotron  period  is  of  the  order  of  twice  the  length.  Wien  the 
signal  amplitude  is  small  compared  to  this  characteristic  amplitude,  one  is 
in  the  linear  regime.  Then  the  relation  between  energy  loss  and  induced 
energy  spread  is  governed  by  the  Madey  Theorem^  (Appendix  A) 

<rf  -  (Yf '  V2  > 

where  <  >  denotes  average  over  entry  phase,  and  implies  that  the  induced 
energy  spread  becomes  much  larger  than  the  mean  energy  loss  as  the  optical 
amplitude  becomes  small.  These  results  have  important  consequences  for  the 
operation  of  an  FEL  as  an  oscillator  in  a  storage  ring.  Renieri’^  Deacon, 
et  al*®  and  Elias,  et  al^  have  shown  that  if  one  depends  upon  the 
synchrotron  radiation  to  damp  the  growth  of  velocity  spread  then  the 


maximum  signal  power  which  can  be  extracted  from  the  electrons  is  of  the 
order  of  the  product  of  the  synchrotron  radiation  rate  with  the  maximum 
value  of  the  fractional  energy  spread  that  will  circulate  in  the  ring, 
hence  of  the  order  of  IZ  of  the  synchrotron  radiation  power. 


4.0 


IMPROVED  ENERGY  EXTRACTION  BY  MEANS  OF  ADIABATIC  DECREASE  OF  THE 


RESONANT  ENERGY 


We  have  seen  in  Section  3  that  the  electron  energy  which  can  be 

extracted  in  passing  through  a  fixed  parameter  wiggler  designed  to  operate 

at  some  specified  z  independent  optical  power  level  is  limited  to  some 

fraction  (~  .5)  of  the  bucket  height  given  by  Eq.  (2.56),  so  that 

(<6y>/Y  )  "  la  a  /u  .  On  the  other  hand,  if  the  wiggler  parameters  are 

r  max  V  w  s 

varied  so  as  to  gradually  reduce  y^  as  the  bucket  moves  through  the 
wiggler,  electrons  which  remain  trapped  in  the  bucket  give  up  an  energy 
approximately  equal  to  the  reduction  in  .  In  this  section  we  discuss 

in  an  exploratory  way  some  characteristic  features  of  this  mode  of 
operation. 


4. 1  Operation  at  Constant  Signal  Amplitude 

A  laser  oscillator,  fitted  with  an  optical  cavity  formed  of  highly 
reflecting  mirrors,  operates  at  steady  state  so  that  the  amplification  of 
the  signal  field  a8  in  passing  through  the  wiggler  just  compensates  for 
the  reflection  loss.  Accordingly,  it  is  a  reasonable  approximation  to 
assume  that  the  electrons  see  a  constant  a8  .  On  the  other  hand,  as  we 
shall  find  that  very  large  a8  Is  desirable,  we  nay  suppose  that  only  a 
small  portion  of  the  resonator  is  occuppied  with  coincident  electron  beam 
and  laser  "micropulses".  On  account  of  Eq.  (2.39),  it  is  also  a  good 
approximation  to  neglect  4k 
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It  is  convenient  for  discussion  and  probably  desirable  as  a  design 
characteristic  to  choose  so  that  ^  is  constant.  Then  Eq.  (2.41) 

immediately  provides  the  relation 
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where 
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The  assumption  of  constant  ♦  together  with  Eqs.  (2.40)  and 
(2.41)  provides  a  useful  constraint  on  the  wiggler  design.  One  sees,  for 
example,  that  if  one  chooses  to  keep  ky  constant,  then  Eqs.  (2.40)  and 
(2.41)  are  satisfied  by  decreasing  aw  linearly.  For  this  choice  the 
bucket  (more  accurately,  the  bucket  minimum)  moves  with  constant  velocity; 
the  electron  energy  loss  comes  entirely  from  the  reduction  in  \i  ,  the 
effective  electron  mass.  Alternatively,  if  one  keeps  aw  constant  then 
the  wiggler  wavelength  decreases  linearly  with  z  .  The  effective  mass 
remains  constant,  but  the  bucket  slows  down.  An  intermediate  case  is 
provided  by  keeping  the  wiggler  field,  bw  -  kySy  ,  constant.  Eq .  (2.41) 
is  again  readily  integrable. 
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Eq.  (4.1)  Implies  the  Inequality 
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Eq.  (4.2)  represents  a  constraint  on  the  length  and/or  optical  power 
density  which  is  required  to  obtain  a  desired  change  in  Yf  .  It  is 
probably  desirable,  therefore,  to  design  the  wiggler  magnet  so  as  to 
minimize  the  right  hand  side  subject  to  a  specified  value^O  of  k^L)  and 
bucket  efficiency,  which  we  define  by  ■  (^  -  Yf)  /  Y±  .  The  required 
relations  are  easily  established  for  the  constant  ky  and  constant  aw 
cases.  One  obtains  the  following  simple  results: 

Constant  k,.,  case: 
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It  is  apparent  from  these  results  that  the  constraint  becomes  severe  as 
approaches  unity,  and  that  it  is  less  severe  in  the  constant 
case.  The  constant  b„  case  turns  out  to  be  intermediate  between  the  two 
(see  Fig.  4.7).  Much  of  the  discussion  for  the  remaining  part  of  this 
section  will  for  definiteness  be  with  reference  to  the  constant 
case.  The  above  discussion  suggests  that  it  is  likely  to  be  one  of  the 
better  methods  of  varying  .  It  also  may  offer  hardware  advantages. 


For  the  system  described  by  Eqs.  (4.3)  and  (4.4)  we  have  from 

Eq.  (4.1) 


1  2  -  rv 

k  La  -  —  a  - —  esc  (4.8) 

w  8  2  b  1  -  a  r 

In  the  absence  of  prebunching,  electrons  will  enter  the  wiggler  at 
arbitrary  initial  phase  ^(0)  .  For  electrons  of  energy  y^  ,  only  those 
for  which 


<  4»(0)  <  w  -  *  =  *2 
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will  be  trapped  in  the  bucket,  where  ^  is  defined  by  Eq.  (2.55).  The 

%  ~ 

fraction  trapped  in  the  bucket  is  therefore  f  =  — - -  .  For  electrons 

r  b  2 n 

with  initial  energy  different  from  Y^  ,  the  trapped  fraction  will  be 

smaller  and  no  electrons  are  trapped  if  the  deviation  exceeds  6y 

max 

Eq.  (2.56).  The  product  f^ri^  can  be  thought  of  as  a  sort  of  idealized 
electronic  efficiency.  Fig.  4.1  provides  a  plot  of  f^  as  a  function  of 
.  It  is  apparent  that  increases  of  idealized  efficiency  either  by 
increasing  or  increasing  f^  drives  the  needed  magnitude  of  Las  to 

larger  values. 

Although  we  are  not  carrying  out  any  serious  study  of  the  effects 
of  transverse  variation  in  this  paper,  we  do  wish  to  take  account  of  some 
obvious  constraints  which  are  thereby  imposed  upon  operating  parameters. 

To  control  the  effects  of  transverse  variation  in  aw  we  write 


k  r  o  (4 .10a) 

we  3  w 

The  quantity  o  is  a  parameter,  which  we  shall  typically  set  equal  to 
w 

one,  which  represents  our  ignorance  of  what  it  should  really  be.  It  may  be 
dependent  upon  other  operating  parameters.  In  addition,  we  assume  the 
oscillator  to  be  provided  with  an  optical  resonator  formed  by  a  pair  of 
mirrors  and  producing  a  Gaussian  beam  centered  in  the  wiggler.  This  causes 
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FIGURE  4.1 


fb,  AS  A  FUNCTION 


an  amplitude  variation  both  transverse  to  and  along  the  optical  beam  of  the 
form 


|as(r,z) | 
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where  f  «  -i  —  ,  and  d  ■  z  -  L/2  •  To  control  transverse  and 

2  c  s 

longitudinal  variation  in  ag  we  write 


and 
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where  o  ,  o  play  a  role  analogous  to  a  .  Combining  these  relations 
e  L  w 

we  obtain 
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which  tells  us  that,  apart  from  the  a  factors,  the  left  hand  side  of 
Eq.  (4.8)  is  proportional  to  the  square  root  of  the  peak  power  circulating 
in  the  optical  resonator.  Returning  to  dimensional  variables  we  obtain  for 
this  power 
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where 


and 
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Equation  (4.11)  provides  us  with  some  working  relations  which  we  find  to  be 
useful  for  further  discussion.  It  suggests  that  there  is  a  minimum 
circulating  power  required  to  reduce  the  Y  of  any  electron  by 
6y  -  Y^  •  We  call  this  the  threshold  power  PT  with  the  working 
definition 

p  _  I  »2c 

*T  2  2 

e 

We  also  note  that  for  fixed  op  and  idealized  efficiency  «  nbfb  ’ 
is  possible  to  minimize  the  right  hand  side  of  Eq.  (4.8)  with  respect  to 
.  We  define  Pffl  as  the  result  of  this  minimization,  with  op  set 
equal  to  one.  This  procedure  provides  us  with  a  useful  estimate  of 
circulating  power  required  for  a  given  idealized  efficiency  and  a 
reasonable  basis  for  selection  of  and  ♦  .  Fig.  4.2a  provides  a  plot 

of  PT  and  Pm  as  functions  of  and  respectively.  Fig.  4.2b 

provides  a  plot  of  the  ^  and  functions  of  which  emerge  from 

the  minimization. 
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Computer  simulation,  to  be  described  later,  yields  electronic 
efficiencies  which  are  of  the  order  of  80  to  90  percent  of  .  There  are 
two  reasons  for  departure  from  the  ideal.  First,  electrons  which  are 
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FIGURE  4.2a 

THRESHHOLD  POWER  FOR  TRAPPING.  PT  ,  AS  A  FUNCTION  OF  BUCKET 

EFFICIENCY  wb  . 

MINIMUM  REQUIRED  POWER,  Pm  ,  AS  A  FUNCTION  OF  IDEAL 

EFFICIENCY  i,  . 


initially  trapped  do  not  necessarily  remain  trapped  for  the  entire  transit 
through  the  wiggler.  An  electron  which  escapes  from  a  bucket  at  some 


z  “  2e  experiences  a  fractional  energy  loss  of 


6Y/Y 


Y(0)  -  Y(zg) 

- TUT) -  <  \  * 


Second,  the  electrons  which  are  not  initially  trapped,  i.e.  those  for  which 

the  initial  phase  is  outside  the  4^  <  f  <  range,  all  tend  to  gain 

energy.  The  average  gain  of  these  electrons  is  of  the  order  of 

i  ( 6y)  ;  both  the  fact  that  there  is  an  increase,  and  its  order  of 

2.  max 

magnitude  can  be  inferred  from  Fig.  2.2. 


Detrapping  appears  to  be  the  more  important  of  the  two  effects. 

Its  occurrence  can  be  readily  related  to  the  adiabatic  invariant  defined  in 

Eq.  (2.49).  For  an  electron  which  enters  at  6y  *  0,  i|<  ■  \j/  we  have 

o 

from  Eqs •  (2.44),  (2.46),  (2.47),  (2.48)  and  (2.49) 
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One  of  the  limits  is  given  by  ^  and  the  other  is  related  to  it  by  (see  Fig. 2. 3) 


cos  +  4\  sin  'l'  “  cos  sin  +  -  cos  t  ^  sin  ^ 

llr  22r  oor 


(4.15) 


In  general  one  has 


^1  c  'J'i  <  <^r  ,  <  <P2  <  <P2  -  it  -  Vf 


so  that  which  of  the  two  is  equal  to  i|>q  depends  upon  the  range  in  which 
'I*  is  located.  In  addition,  at  any  z  there  is  a  maximum  possible  value 
for  J  which  is  given  by 
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The  function,  a(i^  ,  it  -  <|>  )  =  a(i{)  )  is  plotted  as  a  function  of  ip  in 
r  r  r  r 

Fig.  2.6.  The  value  (ze)  of  z  at  which  a  particle  escapes  from  the 
bucket  is  determined  by 
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For  the  oscillator  problem  which  we  have  been  studying,  with  constant 

a  ,  ,  k  and  6k  «  0  we  have 
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and  hence,  from  Eq.  (4.4)  and  (4.5) 
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One  easily  sees  that  for  a  given  ,  there  is  a  range  of  \|;  about  <J»r 

such  that  Eq.  (4.18)  yields  values  of  z£  >  L  .  These  correspond,  of 
course,  to  particles  which  do  not  detrap.  Fig.  4.3  provides  a  family  of 
2  2 

plots  of  fci  (i|>r,  ^Q)/ct  (<|^)J  as  a  function  of  i|i  for  various  values  of 

♦  ,  and  Figs.  4.4a  and  4.4b  present  ze  as  a  function  of  \|;  with 

and  i|»  as  parameters.  Incidentally,  the  adiabatic  theory  of  detrapping 

typically  implies  an  upper  bound  of  electronic  efficiency  for  some 

designs.  For  the  constant  k  ,  'p  ,  a  wiggler  the  bound  is  found  to  be 

w  r  s 

.672.  If  a^  rather  than  ky  is  held  constant,  the  bound  is  reduced  to 
.544. 


A  numerical  simulation  of  the  constant  k  ip  a  case  has 

v  r  s 

been  carried  out,  based  upon  Eq.  (2.25).  The  solution  is  used  to  compute 
(L)  from  Eq.  (2.24)  as  a  function  of  for  various  values  of  and 

♦  •  The  results  of  some  of  these  computations  are  shown  in  in  Fig.  4.5 
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THE  DETRAPPING  FUNCTION  Rj  lip0.<pr)  =  a2(ipr  ip0)/a2(ipf) 
PLOTTED  AS  A  FUNCTION  OF  INITIAL  PHASE  ip  Q  FOR  VARIOUS 
VALUES  OF  RESONANT  PHASE  if,  . 


FIGURE  4.5 


ENERGY  SPECTRUM  OF  EMERGING  ELECTRONS  AS  A  FUNCTION  OF  ENTRY 
PHASE.  ALL  ELECTRONS  ENTER  WITH  ENERGY  EQUAL  TOTr(0). 

THE  THREE  CURVES  REFER  TO  MINIMUM  POWER  CONFIGURATIONS  OF 
THE  SORT  ILLUSTRATED  IN  FIGURE  4 2.  THE  PARAMETERS  OF  THE 

THREE  CURVES  ARE: 
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31.9° 

.215 

.181 
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.66 

21.9° 

.400 

.314 

c 

.823 

11.4° 

.600 

.448 

THE  VERTICAL  LINES  ABOVE  THE  CURVES  DENOTE  POSITIONS  OF  THE 
TRAPPING  BOUNDARIES  AND  ■ 


for  some  of  the  minimum  power  configurations  of  Fig.  4.2.  The  flat 
portions  at  YfU/Y^  ■  1  -  corresponds  to  the  range  of  i|i  over  which 
there  is  no  detrapping. 

The  fact  that  the  rapidly  changing  portions  are  not  vertical  is  an 

indication  that  detrapping  is  taking  place.  The  peculiar  jog  at  the  right 

hand  end  just  before  ^  reaches  is  due  to  particles  which  detrap 

after  the  first  synchrotron  oscillation.  Note  that  untrapped  particles 

gain  energy.  All  of  the  curves  shown  are  for  Y^  *  Yf  ,  with  no  energy 

spread  included.  They  should  be  representative  for  a  spread  small  compared 

to  6y  .  The  numerical  simulation  also  provides  a  value  for  the 
max  r 

electronic  efficiency  ne  .  Fig.  4.6  shows  the  extent  to  which  ng 

differs  from  for  the  configurations  of  Fig.  4.2. 

It  is  clear  from  Eq.  (4.16)  that  detrapping  can  be  avoided  for  the 
constant  ag  case  by  designing  the  wiggler  so  that  aw/kw  Is  constant. 

As  in  the  constant  1^  t  constant  ay  ,  and  constant  bw  cases  discussed 
previously  one  can  use  the  condition  of  constant  4^  and  specified  to 

determine  the  z  dependence  of  aw  and  kw  »  an<*  the  required  power.  The 
required  threshold  powers  in  the  constant  k^,  aw  ,  bw  ,  and  a^/k,,  cases 
are  compared  in  Fig.  4.7  as  functions  of  .  Since  the  penalty  due  to 
detrapping  is  rather  small,  one  actually  obtains  a  better  net  efficiency 
from  a  constant  wiggler  than  from  one  which  avoids  detrapping.  Some 

small  improvement  can  probably  be  obtained  over  the  constant  kw  case  by 
taking  detrapping  into  account,  but  we  have  not  pursued  this  question 
further.  This  assessment  applies  to  the  case  in  which  the  incident 

electron  energy  spread  is  small  compared  to  ^max  •  When  they  are 

comparable,  as  is  likely  to  be  the  case  In  practice,  the  penalty  imposed  by 
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FIGURE  4.7 


COMPARISON  OF  THRESHOLD  POWERS  AS  FUNCTIONS  OF  V  ^  FOR 
VARIOUS  CONSTANT  ipr  CONSTANT  a,  WIGGLER  DESIGNS.  THE 
FOUR  CASES  ARE: 

•)  Constant  aw/kw 

b)  Constant  ayy 

c)  Constant  bw 

d)  Constant  1^ 


detrapping  Increases,  and  Its  avoidance  acquires  more  significance  as  a 
design  consideration.  It  should  be  noted  that  In  the  case  of  an  amplifier 
or  heavily  loaded  oscillator,  where  ag  increases  with  z  ,  detrapping 
tends  to  be  avoided  automatically,  as  discussed  later  In  this  section. 

In  all  of  the  previous  discussion  we  have  imposed  the  requirement 
(largely  to  simplify  the  discussion)  that  ^  be  z  independent.  One  may 
ask  whether  this  choice  leads  to  an  approximately  optimized  wiggler 
design.  From  the  point  of  view  of  detrapping,  it  would  be  useful  to  have 

decrease  with  z  ,  to  compensate  for  the  decrease  in  aw/kw  which 
seems  to  occur  for  the  more  optimum  constant  ^  designs.  This  is 
apparent  from  the  behavior  of  0(1]^)  (see  Equations  (4.13),  (4.14)  and 
(4.18)]  shown  in  Fig.  2.6.  On  the  other  hand,  from  the  point  of  view  of 
initial  capture,  it  would  be  best  to  have  ^  have  a  very  small,  or  even 
vanishing  initial  value.  This  can  be  seen  from  the  plot  of  the  trapping 
fraction,  f^  ,  shown  in  Fig.  4.1.  Thus  to  maximize  trapping  ^  should 
increase  with  z  ,  while  to  avoid  detrapping  ♦  should  decrease  with 

z  .  One  can  in  principle  increase  aw/kw  to  overcome  the  effect  of 
decreasing  t|»  and  thus  avoid  detrapping,  but  experience  with  the  constant 
♦  case  suggests  that  this  will  not  improve  matters  either.  Hence  the 
constant  ,  constant  design  is  probably  not  very  far  from  an 

optimum  one.  Again  in  anticipation  of  the  amplifier  case,  we  note  that  an 
increase  of  ag  with  z  can  eliminate  the  detrapping  associated  with  an 
increase  in  ^  ,  so  that  an  increase  of  with  z  should  have  some 

advantages  for  an  amplifier. 
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We  conclude  this  discussion  by  listing  the  properties  of  a 
"standard"  wiggler  design,  based  upon  the  constant  kw  case.  It  is 


certainly  not  unique  and  probably  not  optimal  either,  but  it  does  provide  a 
convenient  point  of  reference  for  future  discussion.  We  begin  by 
specifying  an  Idealized  efficiency  ,  a  signal  frequency  ojs  ,  and  an 
electron  beam  radius  r_  •  The  associated  bucket  efficiency,  n.  >  and 
resonant  phase,  <|»r  ,  are  determined  from  n  by  Figure  4.2.  For  wiggler 
design  parameters  we  then  propose 
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the  electron  beam  energy 
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is  given  by 
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(4.23) 


The  required  peak  beam  current  is  given  by 


eP„ 

I - B—  (4.24) 

e  mcZYineQ 


where  Q  is  the  "Q"  of  the  optical  resonator,  and  ng  <  ,  the 

electronic  efficiency.  Rewriting  (4.24)  yields 


_ 1 _ 

Q(6wsre/c) 1/2 


<2"V‘ 


(2-2v^) 


1/2 


(4.25) 


where 


3 

17.0  Kilo  Amperes  . 

The  degradations  associated  with  the  variation  of  the  resonator  mode  over 
the  electrons,  as  well  as  the  effect  of  energy  spread  of  the  incident 
electrons  have  been  neglected  in  the  above.  For  assessing  the  effect  of 
the  latter  we  note  the  relation  (using  4.8  and  2.56  together  with  the 
above)  . 
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The  electron  beam  radius  which  can  be  used  is  primarily  determined 
by  the  electron  beam  emittance.  It  is  shown  in  Appendix  C  that  if  one 
assumes  ideal  emittance  scaling,  then  the  beam  radius  must  exceed  a  lower 
limit  proportional  to  the  emittance  constant.  Finally,  we  repeat  that  we 
have  neglected  the  possiblity  that  some  modification  in  the  parametric 
dependencies  may  emerge  if  the  omitted  "o"  factors  are  taken  into  account. 

It  is  instructive  to  illustrate  the  above  expressions  with  a 
reasonable  numerical  example.  Let  us  suppose  we  wish  to  construct  a  1 
micron  oscillator.  A  preliminary  perusal  suggests  that  *  .215  might 
yield  parameters  achievable  with  current  technology.  Considerations 
bearing  upon  the  choice  of  electron  radius  will  be  discussed  in  Appen¬ 
dix  C,  but  for  the  moment  we  simply  specify  re  “  1*24  mm.  With  these 
choices  the  following  numbers  emerge: 


ne 

- 

.18 

“ 

1.75  -  1.18 

z/L 

m 

.43 

Bw 

- 

4.56  -  3.07 

z/L  KG 

*r 

- 

32° 

Ee 

- 

111  Mev 

Pm 

- 

21  GW 

QTe 

- 

1100  Amps. 

X 

W 

m 

2.3  cm 

6W* 

- 

.014 

L 

m 

10.9  m 

As  indicated,  the  beam  current,  and  accordingly  the  optical  power  generated 
(p  /Q,  of  course),  will  depend  upon  the  resonator  Q  which  is  chosen, 
and,  on  account  of  mirror  losses,  the  optical  power  which  emerges  will  be 
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reduced  from  that  generated.  Furthermore,  we  again  emphasize  that  these 
are  peak  values.  If  the  electron  beam  is  provided  by  an  RF  linac,  "peak" 
refers  to  values  which  obtain  during  a  micropulse.  Average  values  during  a 
macropulse  are  reduced  by  a  factor  which  typically  lies  between  25  and  100. 


4. 2  Operation  with  Amplifying  Signal  Amplitude 

In  discussing  the  behavior  of  the  standard  design  the  increase  in 
the  signal  amplitude  as  was  neglected.  This  neglect  is  justified  for 
large  values  of  Q  .  As  one  decreases  Q  we  expect  hg  to  approach 
n  more  closely,  but  otherwise  the  basic  formulas  for  the  operation  of  the 
device  are  still  expected  to  hold  approximately  despite  the  fact  that 
as  is  increasing.  The  decrease  in  Q  is  likely  to  be  associated  with 
increased  output  coupling.  Ihis  fact,  together  with  the  increase  in 

means  that  the  extra  beam  power  input  associated  with  the  increase  in 
le  is  being  more  efficiently  converted  to  signal  power  output. 


These  expectations  are  easily  understood  by  making  reference  to 
Eq.  (2.41),  which  we  rewrite  as 


a  sinty 
8  r 


w  a  dz 
s  w 


(4.27) 


The  right  hand  side  of  this  equation  has  already  been  determined 
by  the  standard  design  to  be  a  constant.  Associated  with  this  design  was 
an  input  value  of  as  and  an  input  value  of  Vr  •  However,  if  a8 
increases,  ♦  simply  decreases.  Reference  to  Eq.  (2.59)  and  Figure  2.6 
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tells  us  that  the  phase  area  of  the  bucket  increases.  Hence  particles 
initially  trapped  remain  trapped  and  apart  from  the  small  increase  in 
energy  of  the  untrapped  particles,  ne  “  ni  “  *  1)116  t0  t)ie  ^act  t'iat 

the  bucket  is  actually  widening,  some  initially  untrapped  particles  may 
even  become  trapped,  and  can  exceed  n  .  The  important  lesson  to  be 
learned  is  the  fact  that  a  more  rapid  than  "planned"  increase  in  as  does 
no  harm  to  the  trapping  and  deceleration  process.  On  the  other  hand  if 
ag  increases  less  rapidly  than  planned,  a  reduction  in  phase  area  occurs 
which  can  lead  to  some  detrapping.  Indeed  if  ag  drops  sufficiently  below 
its  planned  value  to  require  sin ^  >  1  .total  detrapping  takes  place. 

It  should  now  be  apparent  that  a  simple  way  to  design  an  amplifier 
is  simply  to  remove  the  mirrors  from  the  standard  design  and  supply  a  beam 
current  given  by  Eq .  (4.25)  with  1/Q  replaced  by  (G-l)  ,  where  G  is  the 
power  gain.  The  required  input  signal  power  is  given  by  Pm  ♦ 

In  our  discussion  of  high  Q  oscillators  we  established  a  quasi 

optimi  ,ed  connection  between  the  minimum  circulating  power  Pm  and  • 

Clearly  similar  optimization  could  be  desirable  for  an  amplifier.  The 
simple  amplifier  designed  above  makes  no  attempt  to  take  a  programmed 
increase  of  as  into  account  in  specifying  its  yr  profile  and  hence  is 

unlikely  to  be  optimum  in  this  respect.  What  we  would  wish  to  have  is  a 

* 

procedure  for  designing  the  wiggler  so  as  to  maximize  the  efficiency  for  a 
specified  power  output  and  specified  gain. 

As  a  first  example  of  a  step  in  the  desired  direction  we  describe 
a  procedure  for  choosing  as  (and  hence  y  )  so  as  to  maintain  an 
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approximately  constant  •  In  order  to  be  specific  we  consider  the 

constant  case.  We  then  find  from  Eqs.  (2.40)  and  (2.41),  neglecting 
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Note  that  the  coefficient  in  brackets  is  a  constant.  To  relate 
to  aw  we  make  use  of  Eq.  (2.26),  assuming  that 
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~TJo) 
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<Y(0)>  -  <Y(z)>  “constant  x(Yr(0)  -  Yr(z))  .  (4.29a) 


This  assumption  is  based  upon  the  idea  that  the  trapped  electrons  represent 


a  constant  fraction  of  the  electrons,  and  that  they  give  up  energy  at  the 

dYr 

rate  .  It  yields  the  equation 
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where  G  is  again  the  power  gain.  Combining  Eqs.  (4.28)  and  (4.29)  yields 
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The  left  hand  side  of  Eq.  (4.30)  can  be  regarded  as  the  analogue  of  the 
right  hand  side  of  Eq .  (4.2),  that  is,  as  a  quantity  which  in  this  instance 
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should  be  minimized  with  respect  to  a^d)  and  aw(0)  subject  to  the 

constraint  1  -  -  n.  .  This  condition  then  determines  the 

W  (0)  D 

appropriate  values  of  a  (0)  and  a  (L)  and  hence  P(0)  and  P(L)  . 

w  w 

Subsequently  a^(z)  is  determined  from 
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As  an  example,  we  consider  the  case  =  .43  and  for  simplicity  neglect 

the  final  optimization.  Instead  we  obtain  aw(0)  and  aw(L)  from  Eq . 

(4.4),  which  holds  for  the  G»1  case.  For  G**100,  the  value  of  the  left  hand 

side  of  Eq.  (4.30)  is  found  to  be  .203.  This  number  for  the  right  hand 

side  of  Eq.  (4-30)  should  be  compared  to  a  (0)  -  a  (L)  =  1.18  ,  which 

w  w 

would  be  required  by  the  simple  amplifier  design  discussed  previously  in 
which  a^  decreases  linearly.  Thus  the  improved  wiggler  profile  permits  a 
factor  thirty  reduction  in  the  operating  power  level.  The  minimum  of  the 
left  hand  side  of  Eq.(4.30)  is  found  to  be  very  broad  for  G“100,  so  that 
the  value  obtained  by  choosing  aw(0),  a^d)  as  noted  above  is  within 
l/2£  of  the  minimum. 


The  predicted  behavior  of  ag  is  obtained  from  Eq.  (4.29)  by 
making  use  of  a^d)  as  given  by  Eq.  (4.31).  To  complete  this  example  we 
apply  the  standard  procedure  to  the  rest  of  the  design.  Hence  we  choose 
sin  from  Figure  4.2  to  obtain  the  value  .53  and  an  of  .215.  The 
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required  optical  input  power  is  then  640  MW.  Assuming  an  electronic 
efficiency  equal  to  and  a  gain  of  100,  we  obtain  a  required  electron 

beam  power  of  295  GW.  Figure  4.8  shows  the  design  behavior  of  aw  as 
given  by  Eq.  (4.31)  and  as  as  given  by  Eq.  (4.29)  (curve  d).  It  also 
shows  the  computed  values  of  ag  and  electronic  efficiency  as  obtained 
from  numerical  simulation.  The  numerical  simulation  was  based  upon 
Eqs.  (2.25)  and  (2.26),  neglecting  the  6kg  terms,  but  explicitly 
carrying  out  the  averages  required  for  Eq.  (2.26).  Three  cases  are  shown, 
corresponding  to  the  design  beam  power,  a  somewhat  higher  power  and  a  lower 
power.  The  design  beam  power  does  not  quite  yield  the  design  performance. 
First  we  note  that  the  expectation  of  the  ideal  electronic  efficiency 
failed  to  take  account  of  energy  gain  from  the  untrapped  electrons. 

Second,  Eq.  (4.29a),  the  assumption  of  a  constant  ratio  between 
Yr(0)  -  Yf(z)  and  <Y(0)>  -  <Y(z)>,  is  clearly  not  satisfied  at  the  input 
of  the  amplifier.  The  flat  behavior  of  ag  at  the  input  reflects  the 
fact  that  bunching  must  take  place  before  amplification  can  begin.  We  see 
that  increasing  the  beam  power  improves  the  efficiency  as  well  as  providing 
more  than  the  design  gain.  The  reduced  beam  power  example  illustrates  the 
effect  of  a  smaller  increase  in  ag  than  was  anticipated  in  the  wiggler 
design.  Equation  (4.27)  causes  to  constantly  increase,  continually 

detrapping  particles,  until  at  z/L  ~  .8,  <1^  has  reached  — .  All 
particles  are  detrapped  and  amplification  has  ceased. 

Because  ag  is  an  increasing  function  of  z  it  should  be 
possible  to  allow  to  increase  with  z  without  detrapping 

particles.  One  simply  requires  that  the  phase  area  J  ,  given  by 
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Qe(z)/Qc(0) 


Z/L 


FIGURE  4.8 


SIMULATED  PERFORMANCE  OF  A  MATCHED  PROFILE  AMPLIFIER. 


Peak  Electron  Beam  Power 

a)  350  GW 

b)  295  GW 

c)  145  GW 


Electronic  Efficiency 
.214 
.203 
.091 


Curve  d  shows  the  as(z)/as(0)  profile  obtained  using 
Eq.  (4.29a). 
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to  vanish  at  the 


d  Yr 

Eq.  (2.59),  not  decrease.  Indeed,  by  choosing  — 
input  one  can  begin  by  trapping  all  particles,  since  will  initially  be 

zero.  Since  as  will  at  each  z  increase  monotonically  with  increasing 
beam  power,  there  will  always  be  a  critical  beam  power  above  which  all 
particles  remain  trapped,  so  that  ffc  “  1  n£  *  •  Whether  it  is 

better  to  use  the  beam  power  to  increase  f^  to  unity  rather  than  to 
increase  is  a  question  which  we  have  not  explored. 

The  amplifier  length  provides  another  quantity  which  can  be  varied 
as  a  part  of  an  optimization  procedure.  In  our  discussion  up  to  this  point 
we  have  assumed  the  length  to  be  given  by  Eq .  (4.20),  but  there  is, 
however,  no  compelling  reason  to  limit  the  length  in  this  way. It  can  be 
argued  that  so  long  as  the  rate  of  amplification  is  sufficiently  great  to 
maintain  as  constant  on  the  electron  beam  (i.e.,  to  counter  the  effect  of 
diffraction  spreading)  the  device  can  in  principle  be  arbitrarily  lonp.  It 
appears  liklely  that  6kg  can  also  play  an  important  role  in  this 
situation.  We  have  not,  however,  carried  out  a  complete  analysis  of  this 
super-long  amplifier  and  will  not  discuss  it  further  here. 

In  our  discussion  of  the  amplifier  we  consistent  lv  neglc  ted 

fik  .  It  would  of  course  be  a  straightforward  matter  to  includt  it  in 
8 

6ks 

carrying  out  simulations.  It  appears  as  an  order  — —  correction  to 

w 

Eqs.  (2.24)  and  (2.25).  Just  as  a  local  value  of  as^(z)  is  obtained  from 
Eq.  (2.26),  a  local  value  of  6ks  can  be  obtained  from  Eq .  (2.36).  For  the 
case  illustrated  by  Figure  4.8  we  consider  that  <cos  4'/Y>  ■  0  for  the 
first  five  percent  of  the  amplifier,  so  that  6ks  is  in  fact  zero 
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initially.  For  the  region  in  which  afl  increases  approximately  linearly 

it  is  likely  that  -  1  so  that  &8  -  £  "  -(T  and  we 

10  6ks  8 

estimate  <5k  <  —  everywhere.  Thus  — —  <  .003  everywhere,  so  that 

S  L  K 

W 

its  neglect  appears  to  be  justified  so  long  as  transverse  effects  are 
ignored. 


We  conclude  this  section  by  pointing  out  the  fact  that  the  neglect 
of  transverse  effects  appears  to  be  especially  suspect  for  the  design  and 
analysis  of  amplifiers.  In  the  amplifier  illustrated  in  Figure  4.8,  and  in 
the  qualitative  discussion  of  other  designs  that  followed  we  have  depended 
upon  the  predicted  growth  of  a8  to  maintain  the  electron  traps.  We  have 
assumed  in  carrying  out  the  design,  and  in  the  simulations,  that  there  is 
no  transverse  variation  in  the  signal  amplitude.  This  certainly  will  not 
be  the  case  in  actuality.  Suppose,  for  example,  that  the  electron  beam  has 
a  Gaussian  profile.  Then  we  may  write  for  Eq.  (2.37) 


as  “  c  aw  <sin  ^Y>  exP  t“r2/ re  )  (4. 32) 

s 

Where  <sin  <JVy>  should  also  be  assumed  to  be  r  dependent.  We  also  have 

Eq.  (4.27),  which  implies  that  a  sin  ^  is  r  independent.  Suppose  the 

s  r 

amplifier  has  been  designed  so  that  ^  remains  constant  in  z  at  some 
value  of  r  ■  rj  .  For  smaller  values  of  r  we  expect  as  to  increase 
more  rapidly,  decreases  and  trapping  is  unaffected.  For  particles 

outside  r^,  aa  increases  less  rapidly  than  programmed.  Detrapping  can  be 
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expected  to  develop  at  larger  r  and  progressively  work  its  way  into 


r<j  •  This  could  lead  to  a  decrease  of  amplitude  at  r<j  itself  and  hence  a 
progressive  erosion  of  the  edge.  A  similar  erosion  may  be  expected  to 
occur  at  the  trailing  edge  of  the  pulse.  Thus  while  electrons  at  the 
trailing  edge  can  initially  be  trapped  by  the  input  signal,  the  failure  of 
the  signal  to  amplify  at  the  rear  edge  will  lead  to  a  progressive 
detrapping  that  moves  up  the  pulse.  These  considerations  suggest  that  some 
of  the  advantages  of  sophisticated  designs  over  the  "simple"  design  (in 
which  one  depends  upon  the  initial  amplitude  to  maintain  trapping 
throughout)  may  be  substantially  reduced. 

The  situation  is  further  complicated  by  the  fact  that 


5k 


s 


0,2(0)  a; 

lire-  TTFTzT  <cos  */Y>  exp 

s  s 


(A. 33) 


The  presence  of  this  variation  tends  to  focus  light  towards  the  center, 

reinforcing  the  tendencies  discussed  in  the  preceding  paragraph.  In  the 

amplifier  illustrated  in  Figure  A. 8  we  estimated  5k  ~  - — —  -  for 

8  L  +  9z 

z  >  .1  L  •  According  to  Eq.  (A. 33)  this  value  should  be  construed  as 

refering  to  the  center  of  the  beam  so  long  as  the  variation  of  ag  an(j 

<cos  with  r  is  neglected.  Because  6k  decreases  with  z  ,the 

s 

focusing  tendency  is  largest  near  the  beginning  of  the  amplifier.  To 
obtain  a  crude  quantitative  estimate  we  consider  the  problem  to  be  one  of 
propagation  of  rays  in  a  lenselike  medium,  with  index  of  refraction  given 
by  5n  *  c^8^w8  *  According  to  the  paraxial  ray  equation,  we  have 


d 


d6n 

dr 


(4.34) 


where  r(z)  is  the  ray  coordinate.  We  confine  our  attention  to  the 
central  portion  and  write 


d6n 

dr 


H  (6">r. 

r 

e 


(4.35) 


Taking  advantage  of  the  fact  that  rgdnQ/dz  is  small  we  obtain,  for  a  ray 
emitted  at  z  -  Zj  with  zero  slope, 


r(z) 


r  (*j); 


fin  (z) 
o 


cos 


dz/r  + 
e 


fin  'r 
o  e 


(2fino) 


3/2 


(4.36) 


For  the  amplifier  of  Fig.  4.8,  using  the  above  mentioned  estimate  of  6kg 
we  find  for  Zj  *  .1  L  that  focus  takes  place  at  z  «=  .71  L  .  This 
presumably  overestimates  the  effect  because  we  have  neglected  the  fact  that 
the  focusing  effect  further  increases  a_  ,  thus  reducing  6k  .  (The 
factor  <cos  ■♦'/ Y>  also  increases,  but  the  size  of  the  increases  is  limited 
by  the  fact  that  |cos  'j'l  <  1  ,  and  the  increase  in  ag  is  the  dominant 
effect).  Nevertheless  this  estimate  suffices  to  indicate  that  the  effect 
can  be  significant.  We  remark  that  Eq.  (4.33)  does  not  necessarily  imply 
that  6k  is  a  maximum  in  the  center  as  the  increase  in  a_  can  dominate 

S  8 
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other  effects.  The  presence  of  this  factor  in  the  denominator  also 
indicates  that  a  true  self  focusing  instability  does  not  exist.  The  effect 
clearly  decreases  when  the  gain  is  smaller,  60  that  the  net  outcome  could 
simply  be  a  limit  on  useable  gain.  On  the  positive  side,  the  focusing 
effect  can  also  lead  to  some  optical  beam  trapping,  which  could  be  helpful 
for  devices  with  L  >  kre^  •  Clearly  the  complex  of  issues  raised  here 
constitutes  an  important  subject  for  further  research. 

As  a  final  comment  we  note  that  for  an  oscillator, 

J  6k  dz  ~  (4.37) 

o 

This  should  be  compared  to  the  total  phase  shift,  2  tan-1  .5  ,  which  occurs 
in  the  traverse  of  the  wiggler  due  to  diffraction.  As  a  consequence,  low 
Q  oscillators  may  require  some  small  optical  correction  to  take  into 
account  both  the  effect  and  the  fact  that  it  may  vary  during  buildup  of 
oscillation.  Apart  from  this  possibility,  the  6kg  effect  does  not  appear 
to  affect  oscillator  operation. 
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5.0  ADIABATIC  CAPTURE,  DECELERATION,  AND  DECAPTURE 

In  this  section  we  will  treat  the  problem  of  capturing  the 
electrons  into  a  stationary  bucket  with  a  resonant  energy  equal  to  the  mean 
electron  energy,  while  at  the  same  time  minimizing  the  increase  in  the 
phase  area  occupied  by  the  electrons.  After  capture  the  average  phase 
angle  of  the  electrons  must  be  changed  to  a  positive  value,  so  that  the 
electrons  can  be  placed  in  the  center  of  a  moving  bucket  and  subsequently 
decelerated.  For  a  single  pass  device  it  will  be  only  necessary  to  make 
sure  that  the  phase  area  of  the  captured  electrons  is  less  than  the  area  of 
the  decelerating  bucket.  However,  for  a  multiple  pass  system  we  want  to 
minimize  the  Increase  in  phase  area  occupied  by  the  beam.  After  the  beam 
has  been  decelerated  the  average  phase  angle  of  the  electron  is  returned  to 
zero  and  placed  in  a  new  stationary  bucket  with  a  resonant  energy  equal  to 
the  final  mean  electron  energy.  Finally  the  decapture  is  performed  in  such 
a  manner  as  to  achieve  the  small  desired  energy  spread.  The  total  purpose 
of  all  of  these  processes  is  to  obtain  a  minimum  increase  in  the  energy 
spread  while  at  the  same  time  to  reduce  significantly  the  average  electron 
energy. 


To  illustrate  these  processes  we  divide  the  wiggler  into  five 
regions:  region  one,  with  0  <  z  <  z^  ,  is  where  the  adiabatic  capture 

occurs;  region  two,  with  z^  <  z  <  z^  ,  is  where  the  average  phase  angle 
is  increased;  region  three,  with  Z2  <  z  <  z^  ,  is  where  the  deceleration 
occurs;  region  four,  with  z^  <  z  <  z^  ,  is  where  the  average  phase  angle 
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,  is  where  the 


is  decreased  to  zero;  and  region  5,  with  z  <  z  <  z, 

4  5 

decapture  occurs.  The  phase  space  occupied  by  the  electrons  during  various 
stages  of  these  processes  is  shown  as  the  shaded  area  in  Figs.  (5.1a  to 
5 . lh)  . 

Region  1 

In  order  to  capture  all  of  the  electrons  without  increasing  the 

phase  area  occupied  by  the  electrons  it  is  necessary  to  use  a  stationary 

bucket,  i.e.  \|>  =  0  ,  and  to  adiabatically  increase  the  height  of  the 

bucket.  From  Eq .  (2.51)  we  see  that  for  *  0  the  bucket  height  6y 

r  max 

is  given  by 


<5y 

max 


2Yr^ 

J 


1  +  a 


(5.1) 


We  want  to  start  with  5'f  «  0  at  z  *■  0  and  increase  <5y  with  z 

max  max 

while  keeping  t|>  -  0  .  This  may  be  accomplished  by  having  aw  increase 
with  z  while  demanding  that  kw  also  increase  with  z  to  maintain  the 
condition 


k  (z) 
w 


(5.2) 


Since  6y  (0)  *  0  and  y  is  a  constant,  equal  to  y  (0)  .  we  must 

max  r  r  ’ 

have 


FIGURE  5.1  (a)-(h) 

SCHEMATIC  BEHAVIOR  OF  THE  7,ip  PHASE  SPACE  DISTRIBUTION  AS  A  FUNCTION  OF  3  DURING  AN  ADIABATIC  CAPTURE,  DECELERA¬ 
TION,  AND  DECAPTURE  PROCESS:  (a.)  INITIAL  DISTRIBUTION,  (b.)  DURING  CAPTURE,  (c.)  DURING  INCREASE  OF  AVERAGE  ip .  (d.) 
AND  (e.)  THE  DECELERATION  PROCESS,  (f.)  DURING  DECREASE  OF  AVERAGE  tp  .  (g.)  DURING  DECAPTURE,  (h.)  AFTER  DECAPTURE 


(5.3) 


a  (0)  -  0  and  k  (0) 

w  w 


k 

2Y 


We  combine  Eqs.  (5.1)  and  (5.2)  and  obtain 


[«Y  (z)] 

L  max  J 


4kV 


-  [k  (z) 

s  s  L  w 

7 


kw<0)J 


kB(z)ku(0) 


(5.4) 


From  Eq.  (5.4)  it  follows  that  the  maximum  value  of  the  bucket  height 
occurs  for  a  value  kw(z)  =  2kw(0)  corresponding  to  a  value  of  aw(z) 
equal  to  one.  It  is  desirable  to  have  a  large  value  for  the  final  bucket 
height  so  we  assume  that  region  one  of  the  wiggler  is  designed  to  Increase 
both  the  magnetic  field  aw  and  the  wave  number  kw  to  the  final  values: 


k  (z  )  -  2  k  (0)  (5.5) 

w  I  w 

a  (z  )  -  1  (5.6) 

W  i 

which  yields  a  final  bucket  height 


(z.) 
m  1 


(5.7) 


Note  that  the  final  value  of  the  fractional  bucket  height. 


5Vzi> 

v 


is 


r 

dependent  upon  only  the  field  strength  and  the  wavelength  of  the  radi¬ 
ation.  In  order  to  capture  all  of  the  electrons  in  the  bucket  it  is 
necessary  for  the  initial  energy  spread  in  the  beam  before  capture  to  be 
less  than  the  final  bucket  area,  i.e. 


76 


2V  *  ^initial  < 


<5^ 


(5.8) 


While  Eqs .  (5.5)  and  (5.6)  give  the  final  values  for  the  wiggler  wavelength 
and  magnetic  field  necessary  to  achieve  the  maximum  bucket  area  for 
capture,  the  exact  speed  and  method  of  achieving  these  final  values  have 
yet  to  be  determined.  In  order  to  satisfy  the  adiabatic  condition,  the 
change  in  ky  and  Sy  must  be  sufficiently  slow.  The  exact  number  of 
phase  oscillations  that  are  needed  for  the  adiabatic  condition  depends  upon 
with  what  accuracy  one  needs  the  phase  area  to  remain  constant  and  will 
require  computer  studies  for  the  motion  of  a  large  number  of  particles. 
However,  as  a  rule  of  thumb  the  capture  length  must  be  long  compared  to  the 
final  phase  oscillation  period,  i.e. 


z 


1 


(5.9) 


From  Eq .  (5.9)  we  see  that  it  is  desirable  to  have  a  low  resonant  energy 
and  a  hig1.  optical  field  in  order  to  have  a  short  capture  length.  The 
exact  form  of  the  change  in  ky  and  ay  has  only  one  constraint  given  by 
Eq.  5.2,  so  we  are  free  to  chose  the  simple  form: 


kw(z)  •=  kw(0)  U  +  (f^)2aj  .  and 


(5.10) 


aw(z) 


(5.11) 


with  a  a  constant. 
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In  order  to  decelerate  the  electrons,  after  the  adiabatic  capture 


proceed  to  perform  synchrotron  oscillations  about  the  new  value  of  Yf 
(see  Fig.  5.1c),  with  the  maximum  change  in  41  occurring  after  one-quarter 
of  a  synchrotron  oscillation  period.  If,  on  the  other  hand,  the  change  in 
the  value  of  Yf  is  greater  than  the  bucket  height,  the  magnitude  of  the 
increase  of  4*  will  not  be  limited  (i.e.  the  electrons  are  detrapped) . 
Because  we  want  to  prevent  an  increase  in  the  phase  area  occupied  by  the 
bunch  we  need  to  change  the  phase  and  energy  of  all  of  the  electrons  in  the 
bunch  in  a  coherent  manner.  This  can  be  achieved  if  all  of  the  electrons 
in  the  bunch  are  near  enough  to  the  new  bucket  center  to  be  in  the  linear 
region.  In  order  to  be  able  to  accomplish  this  it  is  necessary  that  the 
phase  area  occupied  by  the  bunch  must  be  small  compared  to  the  bucket  area 
and  the  change  in  the  bucket  center  must  be  smaller  than  its  height.  This 
second  condition  yields 


(f/2)Y  <  <$Y 

r  m 


(5.15) 


If  the  restriction  of  Eq .  (5.15)  is  satisfied,  the  average  values 

of  ♦  and  Y  are  equal  to  the  central  values,  and  we  find  that  after  one- 

quarter  of  an  oscillation  4*  reaches  its  maximum  value 


(5.16) 


while  the  average  energy  of  the  bunch  has  changed  by 


Ay  -  -  f/2  Yr(0) 


(5.17) 
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Region  3 

With  the  new  values  of  4>  and  Y  it  is  now  possible  to  use  a 
tapered  wiggler,  as  was  discussed  in  detail  in  Section  4,  to  produce  a 
decelerating  bucket  and  to  decrease  the  energy  further.  The  energy  of  the 
bunch  has  already  been  decreased  during  the  process  of  increasing  +  and 
if  the  desired  energy  extraction  is  sufficiently  small  it  may  be  possible 
to  skip  the  decelerating  region  entirely. 

Region  4 

In  this  region,  with  z ^  <  z  <  z ^  ,  we  must  reduce  the  average 

phase  angle  to  zero  in  preparation  for  the  decapture;  as  an  example  we 
choose  a  program  similar  to  that  used  in  region  2.  That  is  we  choose  ky 
and  aw  such  as  to  produce  a  stationary  bucket  with  a  resonant  energy 
equal  to  Y(z^)  and  allow  the  bunch  to  perform  one-quarter  of  an 
oscillation  so  that  4*  returns  to  zero  (see  Fig.  5. If).  The  new  average 
phase  angle  and  energy  at  z  »  z4  is  given  by 

>Kz, )  =  0  and  (5.18) 

4 

Y(z.)  -  (1  -  f/2)  Y(z  )  (5.19) 

4  3 

Note  that  if  major  decelerating  region  3  is  eliminated  the  total  energy 
change  is 


AY  -  Y(z4)  -  Y (0)  “  -  fYr 


(5.20) 
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which  corresponds  to  jumping  the  stationary  bucket  down  in  energy  by 
-  f/2  and  allowing  the  bunch  to  perform  one-half  of  a  phase 
oscillation. 

Region  5 

In  the  fifth  and  final  region,  with  z^  <  z  <  z,.  ,  the  electrons 

are  debunched  by  allowing  the  magnetic  field  to  decrease  to  zero.  This 

process  is  the  reverse  of  the  capture  process  of  region  1  and  also  must  be 

done  adiabatically  in  order  to  conserve  the  electron  phase  space  density. 

It  follows  that  the  wiggler  period  must  increase  such  that  'P  remains  zero 

and  ■  constant  .  First  the  center  of  the  bucket  must  be  placed  at  a 

value  of  Y^  «  Y(z^)  by  a  discontinuous  change  in  the  wiggler  period  and 

field  such  that  for  z  >  z, 

A 

k  (z)  =  (1  -  f)  k  (z,)  (5.21) 

w  w  A 

a  (z)  «  a  (z  )  (5 .22) 

w  w  A 

Next  we  must  decrease  the  bucket  area  to  zero  to  decapture  the  electrons; 
for  simplicity  we  choose  the  reverse  of  the  form  chosen  for  the  adiabatic 
capture  in  region  1,  i.e.  for  z^  <  z  <  z^ 


and 


k  (z) 
w 


kw(z4> 


7 


a  (z)  - 
w 


a 


(5.23) 


(5.2A) 
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It  is  clear  that  an  arbitrarily  long  wiggler  magnet  could  in  this 
way  decelerate  electrons  while  producing  little  energy  spread.  However,  in 
order  to  determine  whether  or  not  the  scheme  is  a  promising  one,  further 
analysis  is  required.  Especially,  a  more  quantitative  estimate  is  needed 
of  the  number  of  oscillations  required  in  the  adiabatic  trapping  and 
detrapping  steps  of  the  process.  Furthermore,  because  of  the  additional 
length  associated  with  these  steps,  a  more  severe  side  band  instability 
problem  (see  Sec.  7)  may  be  anticipated. 
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PHASE  AREA  DISPLACEMENT 


In  order  that  a  large  fraction  of  the  electrons  be  decelerated  and 
hence  transfer  energy  into  the  optical  field  for  the  operational  inodes  of 
the  PEL  described  In  the  previous  sections.  It  was  necessary  that  the 
energy  spread  of  the  incoming  electron  beam  be  less  than  the  maximum  bucket 
height.  In  this  Section  we  will  discuss  the  method  of  phase  area  displace¬ 
ment^  which  can  allow  all  of  the  electrons  to  be  decelerated  even  when  the 
initial  energy  spread  (or  effective  energy  spread  when  transverse  emittance 
and  magnetic  field  variation  with  beam  size  are  included)  is  considerably 
larger  than  the  bucket  height.  Phase  area  displacement  refers  to  an  opera¬ 
tional  mode  In  which  an  empty  bucket  is  accelerated  through  the  phase  area 
of  the  beam  with  the  result  that  the  phase  area  occupied  by  the  electrons 
is  displaced  downward  in  energy.  Consider  the  case  wliere  the  accelerating 
bucket  starts  with  a  resonant  energy  far  below  the  energy  of  the  electrons 
in  the  beam  and  is  adiabatically  moved  through  the  beam  until  the  final 
resonant  energy  is  far  above  the  electron's  energy,  as  illustrated  in  Fig. 
6.1.  Note  that  for  an  accelerating  bucket  ^<0  .  The  final  mean  energy 

of  the  electrons  is  lowered  by  the  phase  area  of  the  empty  accelerating 
bucket  divided  by  2 n  ,  while  the  final  energy  spread  of  the  beam  is  nearly 
equal  to  the  initial  energy  spread,  i.e.. 


<Y(0)>  -  <Y (L)> 


(6.1) 


Ay 


f 


Ay, 


(6.2) 
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FIGURE  6.1 


POSITION  OF  EMPTY  BUCKET  AND  PHASE 
AREA  OF  ELECTRONS  AT  VARIOUS  POSITIONS 
IN  THE  FEL. 
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where  J  is  the  area  enclosed  by  the  accelerating  bucket  and  is  given  by 


J 


Eq.  (2.59).  From  Fig.  6.1  we  see  that  the  total  change  in  must  be 

much  larger  than  the  sum  of  the  bucket  height  and  the  energy  spread  in  the 
beam,  i.e. 


[y  (l)  -  y  (0)]  ~  y'  l  »  2  ($y)  +  ay  , 

L  r  r  r  max  spread 


(6.3) 


As  long  as  both  the  adiabatic  condition,  which  is  discussed  below, 
and  tht  above  condition  are  met,  the  energy  spread  of  the  beam  is  not 
greatly  increased  and  the  average  energy  loss  of  the  electrons  is  indepen¬ 
dent  of  the  initial  energy  spread.  We  note  that  the  Hamiltonian  for  the 
particle  motion  as  given  by  Eq.  (2.44)  is 


~  (6y)2-  ~8c^~  (cos  4-  +  4»  sin  4»r) 
r  r 


(6.4) 


It  is  convenient  to  use  a  transformation  of  the  independent  variable 


2k 


du 


dz 


(6.5) 


to  define  a  new  Hamiltonian 


H5(V,6y,u) 


(«Y) 


2  -  A  [cos  +  4>  sin  ) 


(6.6) 


with 


A 


at 


a 


s  w 

2  k  c 
w 


(6.7) 
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If  is  a  constant,  less  than  zero  for  an  accelerating  bucket,  then  for 

the  wiggler  design  with  -r—  ■  constant,  the  Hamiltonian  is  a  constant  of 


w 

the  motion;  for  simplicity  we  will  specialize  to  this  case  to  derive  an 
analytical  expression  for  the  rms  energy  spread  for  the  beam  after  passing 
through  the  FEL. 


The  equation  for  the  rate  of  change  in  the  synchronous  energy  in 
terms  of  u  can  be  obtained  from  Eqs.  (2.41)  and  (6.5) 


—  Y  -  -  A  sin  * 

du  r  r 


(6.8) 


The  rate  of  change  for  the  total  energy,  Y  *  (Ay  +  Yf]  ,  follows  directly 


from  Hamilton's  equation 
Eq.  (6.8)  yields 


3H5 
~ 7$ 


which  when  combined  with 


!u  Y  ■  -  A  sln * 


(6.9) 


We  can  use  the  fact  that  H ^  is  a  constant  of  the  motion  along  with 


3H, 


Hamil 


dii  5 

Tiilton's  equation  ^ t0 


(6.10) 


and 


6y  -  +/2  H  +  A  (cos^  +  'isin4'r,]|1 


/2 


(6.11) 
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Thus  the  change  in  energy  is  given  by 


Ay 


A  r _ +  (sin 4')  dij1 _ 

/2  J  F - 

Ja  +  A[cos4'  +  ipsin  ^  ) 


(6.12) 


We  are  considering  an  electron  which  passes  completely  around  the  bucket, 
i.e.,  41  starts  at  a  large  negative  value  i|^,  reflects  at  ^  and 

returns  to  a  large  negative  value  ^  as  shown  in  Fig.  6.2.  We  neglect 
the  quantities  which  oscillate  rapidly  as  ^  and  ^  _  «  ,  neglect 

^  compared  to  ^  and  ^  ,  and  note  that  H  -  -  A  [cos  sin  ) 

to  find 


Ay  ■  -  A/2 


% 


/ 


sin^  dv 


[_n/h  +  A  (cos  +  Vsin  ^  J  \i  H  +  A[cos'i<f  +  ^sin^J 

(6.13) 


If  we  know  the  probability  of  an  electron  having  a  certain  value  of  ^ 
then  we  can  use  Eq.  (6.13)  to  calculate  both  the  average  energy  change  of 
electrons.  Ay  ,  and  the  final  rms  spread. 

The  evaluation  of  Eq.  (6.13)  is  rather  complicated  and  details  are 
given  in  Appendix  B.  A  further  energy  spread  is  expected  from  the  fact 
that  the  bucket  is  of  course  not  moved  up  from  to  +°°  but  has  a 

finite  energy  traverse 

2  A~  -  Yr(L)  -  Yr(0)  .  (6.14) 
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The  additional  energy  spread  due  to  the  end  points  is  uncorrelated  with  the 
spread  due  to  the  spread  calculated  above  and  the  two  spreads  should  be 
added  in  quadrature.  This  effect  also  is  considered  in  Appendix  B  where  it 
is  shown  that  for  4*r<<!  «  the  average  energy  gain  is  given  by 


(6.15) 


and  the  ratio  of  rtns  energy  spread  to  energy  gain  is  approximately  given  by 


(AU^ng.  „  /  sln2^  +  __l!  /H  V  .  (6.16) 

T?  V  r  (64) 2  y  Ay  / 


These  equations  are  in  good  agreement  with  numerical  calculations  for  small 
values  of  i|i  . 

For  a  reasonable  wiggler  it  is  not  feasible  to  allow  a  very  large 
value  of  Ay  or  a  very  small  value  of  ,  and  because  of  the  approxi¬ 
mations  made  in  deriving  the  Hamiltonian  it  is  necessary  to 

numerically  integrate  the  equations  of  motion,  Eqs.  (2.22  and  2.23),  for  a 
large  number  of  electrons  with  various  initial  conditions.  For  the  case 
where  ^  is  held  constant  a  small  fraction  of  the  particles  are  often 
captured  by  the  moving  bucket  and  will  have  their  energy  increased,  thereby 
significantly  contributing  to  the  increase  in  the  energy  spread  of  the 
final  electron  beam.  W.’ther  this  trapping,  which  would  not  occur  for  a 
static  Hamiltonian,  is  due  to  numerical  errors  or  the  small  non-static 
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corrections  to  the  Hamiltonian  is  not  now  completely  clear.  However  it 
would  probably  be  present  due  to  small  errors  in  any  physical  realization 
of  the  wiggler.  The  particles  that  may  become  captured  are  those  which  are 
near  the  unstable  fixed  point  of  the  bucket  when  the  bucket  passes  through 
the  phase  space  of  the  beam.  One  method  that  may  be  used  to  prevent  this 
capture  is  to  increase  the  magnitude  of  the  resonant  phase  angle  Yr  as 
the  bucket  is  accelerated;  this  moves  the  unstable  fixed  point  slightly, 
preventing  the  trapping  of  the  electron.  This  results  in  electrons  with 
larger  energies  being  decelerated  less  than  those  at  the  lower  energies 
producing  an  extra  source  of  energy  spread  in  the  emerging  beam.  A  balance 
between  all  of  these  effects  as  well  as  a  desire  to  minimize  the  length  of 
the  wiggler  requires  computer  simulation  to  arrive  at  a  reasonable  design. 

While  this  method  will  allow  a  much  larger  entrance  bandwidth  and 
produce  a  relatively  small  additional  energy  spread  as  compared  to  a 
constant  wiggler  one  must  pay  the  price  with  a  longer  wiggler  magnet.  It 
also  appears  that  the  small  signal  gain  of  such  a  device  behaves  differ¬ 
ently  from  that  of  the  constant  parameter  FEL.^  At  this  writing  it  does 
not  appear  to  be  out  of  the  question  that  the  phase  area  displacement 
wiggler  has  adequate  charac teristics  to  be  used  in  a  CW  storage  ring, 
although  much  more  simulation  is  required  and  the  start-up  problem  must  be 
addressed . 
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V 


7.0  INSTABILITIES 


7.1  Introduction 

We  have  been  concerned  In  this  paper  primarily  with  free  electron 
lasers  operating  at  very  high  Y  ,  and  our  principal  interest  has  been  the 
mode  of  operation  in  which  the  electrons  are  trapped  In  decelerating 
potential  wells  for  a  single  pass  through  the  wlggler.  In  this  mode  of 
operation  collective  effects  should  be  small*  In  particular  the 

/  2\1/2 

relativistic  longitudinal  plasma  frequency  u>  ■  (  ^ne  j  i8 

PL  V  mY3  / 

pLdt  «  1  for  a  pass  through  the  wlggler.  Thus  we  will 
concentrate  our  attention  on  Instabilities  which  arise  essentially  from 
interactions  of  Individual  electrons  with  the  electromagnetic  field. 

These  are  of  two  kinds.  First  there  are  problems  arising  from  the 
transverse  structure  of  the  beam  which  occur  because  the  effective 
refractive  index  depends  on  the  density  of  trapped  electrons.  A 
qualitative  discussion  has  been  given  in  Section  4;  a  more  complete 
quantitative  treatment  is  in  progress  and  will  be  published  later. 

A  second  type  of  interaction  which  can  occur  is  the  unstable 
generation  of  parasitic  electromagnetic  waves,  l.e.  waves  with  a  different 
frequency  from  that  of  the  signal  which  we  are  trying  to  grow.  Such 
generation  can  be  expected  in  any  situation  other  than  the  simple  linear 


enough  that 


/" 
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regime  operated  at  the  maximum  of  the  signal  gain  curve.  Thus  in  the 
"descending  bucket"  scheme,  a  highly  non-linear  mode  of  operation,  the 
electrons  are  all  effectively  trapped  and  moving  on  the  average  at  a  common 
velocity.  We  might  therefore  worry  that  a  wave  of  somewhat  lower  frequency 
than  the  signal  wave,  corresponding  to  a  ponderomotive  phase  velocity 

v  ■  — —  y  which  is  slightly  smaller  than  the  bucket  velocity,  could  be 
R  w 

unstably  amplified.  In  the  remainder  of  this  section  we  will  show  that 
such  amplification  does  occur. 


7.2  Linear  Theory  of  Raman  Instability 

The  problem  which  we  will  treat  is  that  of  the  stability  of  an 
infinitesimal  electromagnetic  wave  as  it  perturbs  the  equilibrium  of 
electrons  trapped  in  the  ponderomotive  buckets  of  the  signal  wave.  In  a 
single  pass  wiggler  this  equilibrium  is  not  quite  a  static  one.  The 
potential  wells  may  change  slowly  in  shape  and  depth  as  the  wiggler  is 
traversed,  the  electron  oscillation  phase  may  not  be  quite  random  due  to 
the  trapping  process  at  the  entrance  to  the  wiggler,  and  the  wiggler  is  of 
course  of  finite  length.  In  order  to  make  the  stability  problem  tractable 
we  will  idealize  to  a  true  static  equilibrium.  This  would  correspond 
physically  to  a  constant  parameter  wiggler  in  which  the  electrons  are 
uniformly  accelerated  by  a  DC  (or  resonant  RF)  longitudinal  electric  field 
of  such  a  magnitude  as  to  exactly  balance  the  radiative  loss  to  the  signal 
wave.  This  equilibrium  may  be  described  by  the  Hamiltonian  H ^  given  in 
Eq.  (2.44) : 
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ku  +  6ka  .  o 

H4  -  - •  (6y)z  -  ■  -cy  (cos  f  +  ♦  sin  t(»r)  (7.1) 

where  we  recall  that  Is  the  resonant  y  for  particles  in  the  well, 

k.,  Is  the  wlggler  wave  number,  6k  the  modification  of  the  signal  wave 

number  which  may  be  lumped  Into  kw  and  will  be  subsequently  Ignored, 

aw  and  as  the  non-dlmenslonal  vector  potentials  of  the  wlggler  and 

z 

signal  fields,  and  the  phase  =  f  (k  +  k  )  dz  -  u  t  .  We  recall  that 

o  w 

the  canonical  variables  in  Eq.  (7.1)  are  6y  and  ^  •  The  Independent 
variable  Is  z  . 


For  the  static  situation  which  we  are  considering  all  quantities 
except  6y  and  $  are  constants  and  from  Eq.  (2.41)  we  see  that  the  D.C. 
electric  field  necessary  to  maintain  the  static  equilibrium  Is  given  by 


awas“s 

V 


sin  ^ 


(7.2) 


We  note  further,  from  Eq.  (2.41),  that  the  rate  at  which  a  given 
electron  is  radiating  into  the  signal  wave  depends  on  Its  phase  and  is 
given  by:  (in  units  of  mc^) 


•  asaw“s 

E  «  — - -  sin  . 

a  T_ 


(7.3) 


In  order  to  eliminate  unnecessary  constants  It  Is  convenient  to 


define  s-^  z^  ,  and  p  -  /  ^w*1  6y  .  In  these  units  the 

cy2  F  ./  UJT_TT 


8  8  W 


Hamiltonian  becomes  simply 
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(7.4) 


2 

H5  -  | - (cos  <|>  +  <|>  sin  ^r) 

where  p  and  i|<  are  the  new  canonical  momentum  and  coordinate  effectively 
and  U  Is  the  Independent  variable. 


This  Is  of  course  the  Hamiltonian  for  the  pendulum  equation  which 

we  have  discussed  in  Sect.  2  and  3  with  a  succession  of  potential  wells  as 

shown  in  Fig.  2.1.  Electrons  trapped  In  the  potential  well  oscillate  with 

frequency  fl  (Hj)  •  This  frequency  ranges  from  «  «  0  for  particles  at 

the  top  of  the  potential  well  to  fl  -  Jcos  ^  for  particles  at  the 

bottom.  Thus  Au  ■  2*  corresponds  to  a  synchrotron  period,  while  a 

—1  /2 

typical  wlggler  length  would  be  of  order  U  (a  )  for  an  efficient 

single  pass  design.  To  specify  the  equilibrium  we  must  give  the  distri¬ 
bution  of  electrons  f(p,<W  in  the  potential  well.  In  steady  state  we 
expect  f  to  be  Independent  of  the  phase  of  oscillation  in  the  potential. 
It  is  therefore  convenient  to  Introduce  action  angle  variables: 


J(H5) 


£  p  d<J»  ■  ji  di>' 


j2(h5  +  cos  V  +  V  sin 


(7.5) 


and 


<J>( 


„r  . 

4  J2  (H5  +  cos  ♦'  + 


sin  ^ ) 


with 


“("s’  ■  2'[*  ^] 


-I 


(7.6) 
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The  symbol  j)  means  a  complete  cycle  of  oscillation  In  the 
potential  well.  From  Eq.  (7.6)  we  see  that  $  varies  between  0  and  2*  in 
an  oscillation  period.  The  area  in  phase  space  is  given  by  dJ  d$  and 
the  equilibrium  distribution  function  is  independent  of  ♦  .  Thus  we 
specify  the  equilibrium  by  giving: 


f  -  f  (J)  .  (7.7) 

o  o 

The  equations  of  motion  of  a  particle  are 


and 


(7.8) 


(7.9) 


Note  that  because  of  the  well-known  adiabatic  properties  of  J,  we  might 
expect  our  results  to  be  applicable  also  to  slowly  varying  wigglers.  For 
convenience  we  normalize 


r 


f  (J)  dJ 
o 


1 


The  rate  at  which  energy  is  given  to  the  signal  wave  may  be  seen 
from  Eq.  (7.3)  to  be: 


asawus 

cY_ 


1 

7W 


f  sin  dJ  d 4> 


(7.10) 
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For  the  equilibrium  distribution  this  is 


dEs 

3z“ 


asawus 


TT 


sin 


(7.11) 


since  <sin  i|fc  ■  sin  i|>r  ,  (the  bracket  Indicating  a  phase  average  i.e., 
<A>  ■  -j^Ad^)  .  Finally  we  note  that  the  phase  space  above  the 

potential  well  must  be  empty,  i.e.,  the  phase  space  density  of  untrapped 
electrons  must  be  much  smaller  than  that  for  trapped  electrons  which  have 
been  held  together  by  the  ponderomotive  wave.  Moreover  we  expect  the 
distribution  to  be  continuous,  since  various  diffusive  processes  which  we 
will  discuss  later  are  operative.  Hence  it  is  reasonable  to  impose  the 
condition 


f  (J  ) 
o  max 


0 


where 


max  corresponds  to  the  top  of  the  well. 


(7.12) 


Next  we  consider  the  introduction  of  a  small  perturbing 
electromagnetic  wave  defined  by  a  helical  vector  potential  analogous  to 
that  given  for  the  signal  wave  in  Eq.  (2.8) 

\  “  £As  [x  cos  (kRz  -  u>Rt)  -  y  sin  (kRz  -  u>Rt)  J  .  (7.13) 

Here  £  specifies  the  relative  strength  of  the  perturbing  wave.  The 
Hamiltonian  H  of  Eq.  (2.9)  then  becomes 
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2  l  2  2 

H'  «  me  )  M  +  t  +  2a  a  cos  y  +  2  e  a  a  cos 
1  )  w  s  v  s  R 


.1/2 


where  the  phase  o£  the  perturbing  ponder omotive  wave  is  given  by 


*R  "  (kw  +  V  Z  “  V  "  kv 


("?) 


(1)D 

2  +-I?  * 

s 


(7.14) 


The  development  of  Section  2  is  easily  modified  to  Include  both 

radiation  fields,  A  and  A  ,  to  yield  the  new  Hamiltonian: 

IT  8 


H5  -  -  (cos  ¥  +  *  8 in  vr)  -  £  COS  VR  , 


(7.15) 


with,  in  these  units,  y  -  v  +  AU  .  The  canonical  variables  of  course 

K. 

remain  p  and  V  *  We  will  retain  the  angle  action  variables  defined  by 
Eqs.  (7.5)  and  (7.6)  although  the  evolution  of  p  and  y  are  now  given  by 
Hamilton's  equations  for  (7.15). 


WR 

a. 

8 


In  this  expression  we  have  neglected  -  i)  y  since 

s 

1  «  1  and  y  remains  bounded  for  trapped  particles.  Further, 


X 


.  1 
'  X 


I1  +  a2 

_asV 


i 


(7.16) 


(We  note  parenthetically  that  for  a  •  1,  X  would  be  nearly  constant  even 

w 

for  a  variable  wiggler  oscillator.  For  an  amplifier,  in  which  as  is  a 
strong  function  of  x,  X  may  vary  appreciably.) 
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Similarly  following  the  developments  in  Section  2  we  find  that  the 
energy  given  to  the  perturbing  wave  is  analogous  to  that  described  by  Eq. 
(7.10),  i.e. 


dER 

dz 


f  sin 


♦R  dJ  d* 


(7.17) 


Since  f  is  independent  of  U  the  energy  transfer  given  by  Eq. 
o' 

(7.17)  will  be  oscillatory  in  U,  i.e.,  vanish  on  average,  for  the 
equilibrium  distribution  f  -  fQ  •  Hence  it  is  necessary  to  proceed  to  the 
next  order  and  determine  fj  f  the  distribution  function  to  first  order  in 
e  ,  the  amplitude  of  the  sideband. 


Noting  that  the  phase  space  element  2*  dpd'l>  2  dJd<j>  ,  we  can  write 
the  Liouville  equation  for  the  distribution  function: 

3f  d$  3f  dJ  3f 

■3Tr+^!I'3?  +  dU'Jj"0 

3f0  3fo 

Recalling  that  -  0  the  linearized  Liouville  equation  is: 


3fi  +  n3fi.dj|  3fo 
-Sir  +  u  +  dul  Tj~ 


0 


(7.18) 


It  is  thus  necessary  to  find  the  rate  of  change  of  J  caused  by 
the  perturbation. 


From  Eqs.  (7.4),  (7.5),  (7.6)  and  (7.15),  we  have: 


dj  2.  iH5  2,  [  w5  w5  ,  ",  1 

du  '  TT  3iT  '  TT  L--3T  “ST"  TF  “STJ 

2  it  Sili 

*  -  Jr  pE  s^n^R  “  "  2*e  T$  sin*R 

-  2*e  4$  cos  *R  . 


We  may  now  substitute  into  Eq.  (7.18) 


2ne  3fO  * 

IT  cos*R  TT  +  f  1 


to  see  that  satisfies: 


3f,  3fi  ,  3f„ 

lu-  +  +  e  ~ir  8in  » ♦)  +  ^ )  -gj-  *  0 


(7.19) 


The  first  term  In  the  above  expression  for  f^  will  give  no 

dER 

non-oscillatory  contribution  to  as  may  be  seen  from  Eq.  (7.17)  since 

cos  sinijL  oscillates  with  frequency  2*U  .  Hence,  only  f 

R  R  *■ 

contributes  to  energy  transfer.  We  solve  Eq.  (7.19)  by  first  substituting 


fj  *  f+  (.!>♦)  e  1  U  +  c.c. 


Further,  writing  f  ■  e  x  It  i»  easy  to  solve  for 
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Tl 


f 


+ 


(r  +4*') 


(7.20) 


The  lower  limit  has  been  specified  by  the  following  argument, 
familiar  in  plasma  physics.  We  must  suppose  the  perturbation  to  be  a 
growing  wave  i.e.  X  has  a  small  positive  imaginary  part.  This  makes  the 
integral  convergent.  Moreover  we  see  that  f+  ( 4>  +  2n)  *  f  (<t>)  as  is 
required  for  a  single-valued  function.  Hence  Eq.  (7.20)  is  the  unique 
solution  of  Eq.  (7.19).  We  may  now  find  the  energy  transfer  rate  from  Eq. 
(7.17). 


dER 

Sz~ 


d*sl"  *r 

r 


-i  Xu 

e  + 


c  .c. 


We  are  in  fact  concerned  with  that  part  of  the  radiated  energy  which  is 
non-oscillatory  in  U  ,  i.e. 


(7.21) 


We  may  simplify  the 
♦  is  a  periodic  function  of 
periods  differ  by  factors  of 


integral  in  the  following  way.  Note  that 
$  so  that  the  integrals  for  successive 
e"lX2,,/n  .  Thus 
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We  may  exchange  labels  <t>  and  <j>"  In  the  second  tern  to  obtain 


2*  4> 


2/  d<t>  f  d<t>' 

n«0  o  o 


2Trm^  -l(r  +4*')  i(4*  +  *) 


e  e 


2iti(n+l)4  “i  +  4^)  +  *') 

4-  e  0  ’  0  ' 


Finally  1  must  be  added  to  Its  complex  conjugate  to  yield 


■E 


F  /■  ,  2*i« 4 ( -1  *'-[♦+ 4  -if*  - v  +4(*-*' 

f  d<t>  I  d<t>  e  le  L  +  e 


n--°°  o  o 


“  2TTin  “  x  ® 

However  /  4  e  -  /  ,  -  m]  -  W6(X  -  mft)  where  m  is  any 

n*-“  m— 1 00  m«-°° 
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integer.  We  may  again  exchange  4>  and  4>'  in  the  second  term  to  obtain 
finally , 


I-E  2  .  (7.22, 

m=-“  o 

This  result  may  now  be  substituted  into  Eq.  (7.21).  However  it  is  useful 
first  to  define  the  relative  gain  of  the  unstable  wave  to  the  signal 
gain.  Thus, 


G(A) 


_l_lVl_^s 

2  2  dz  7  2  dz  * 
a  a 

s  s 


From  Eqs .  (7.11,  7.22  and  7.23)  we  have 


(7.23) 


Equation  (7.24)  is  our  principal  result.  We  note  that  in  an 
amplifier  G  directly  gives  us  the  number  of  e-foldings  of  an  unstable 
mode  relative  to  the  signal  gain.  In  general  it  will  turn  out  that  G  is 
only  slightly  larger  than  1  and  therefore  the  buildup  from  noise  is  not 
significant.  Further  we  see  from  Eq.  (7.16)  that  since  in  an  amplifier 
ag  varies  strongly  with  U  ,  the  resonant  condition  \  -  is  only  met 


for  a  short  time.  We  conclude  that  the  instability  is  not  important  for  an 
amplifier. 


For  an  oscillator,  however,  where  the  optical  waves  are  reflected 
many  times  we  must  require  G(X)  <  1  if  the  external  circuit  is  non- 
dispersive,  i.e.,  if  the  external  fractional  attenuation  of  the  signal  and 
the  sideband  are  equal.  For  G  >  1  the  sideband  would  grow  when  the  sig¬ 
nal  is  in  steady  state,  and  as  we  discuss  in  the  following  section,  would 
presumably  destroy  the  particle  trapping.  Hence  an  effective  requirement 
for  oscillator  stability  is  G(X)  <  1  for  all  X  . 

The  form  of  Eq.  (7.24)  immediately  suggests  the  physical  inter¬ 
pretation  of  our  results.  Recall  that  the  action  J  is  analogous  to  the 
quantum  level  of  the  electron  in  the  potential  well.  Thus  the  instability 
consists  of  a  stimulated  Raman  scattering  in  which  the  signal  wave  uig 
decays  into  the  sideband  plus  m  "quanta"  of  particle  oscillation. 

Since  both  this  process  and  its  inverse  are  possible,  the  net  energy 

3f 

transfer  depends  on  the  derivative  of  the  distribution  function  .  If 

3f 

tj-j—  <  0  ,  the  usual  case,  it  is  waves  with  X  >  0  ,  i.e.,  the  lower  side¬ 

band,  which  grow.  When  many  waves  of  different  frequencies  X  are  excited 
we  may  expect  the  particles  to  diffuse  in  the  potential  well  leading 
eventually  to  detrapping  as  will  be  discussed  later. 

Equation  (7.24)  is  of  course  a  rather  complicated  one,  with  the 
gain  depending  in  detail  on  the  distribution  fucntion.  We  have  not  quite 
been  able  to  prove  that  it  is  impossible  to  make  G(X)  <  1  for  all  X  . 
However  the  following  sum  rule  indicates  strongly  that  this  is  the  case. 
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Consider 


f  V  C  3fo  2  I  f  i('*'+  2  r  i-l 

1  GX  dX  -  -  2^ JdJ  ~  n>  0  I  d*  e  (4  sin  J 

o  m»0  o 


(7.25) 


Now 


i  ( 4*+‘*n  4*) 


2 


o 


E  /«/  «*  3  {£  *1<f  *  ° 

BF-*  O  O 


♦  ') 


/ 


3v|i 

3 


d<i> 


J 

"  n 


In  the  last  line  we  have  used  definitions  (7.5)  and  (7.6). 

Substituting  this  result  into  Eq.  (7.25)  and  integrating  by  parts  on  J 

while  using  the  fact  that  f  (j  )  ■  0  ,  we  find  that 

o  max 


IT 

4  sin  ♦ 

r 


(7.26) 


We  nay  compare  this  to  r  ,  the  value  which  we  would  have  for  a 

1  /2 

marginal  ce.se  where  G  =  1  from  X-0  to  X-fl  «  (cos  ^  )  and 

max  r 

2ero  for  X  >  ft  .We  see  that 
max 


104 


This  result  is  not  quite  conclusive  proof  of  instability  since  G  is  non¬ 


zero  even  for  X  >  ft  due  to  harmonic  emission.  However  in  practice 

max 

these  growth  rates  appear  to  be  small,  and  are  not  enough  to  offset  the 
numerical  factor  in  Eq.  (7.27).  We  would  expect  therefore  that  values  of 
G  somewhat  greater  than  unity  would  typically  be  found.  Ihis  is  indeed 
true  for  all  distributions  we  have  looked  at. 

There  is  a  particular  distribution  function  f  ,  of  some  physical 
interest,  for  which  Eq.  (7.24)  simplifies  somewhat.  This  distribution  is 
given  by 


for  J  <  J 


max 


max 


and 


f  «  0  for  J  >  J  .  (7.28) 

o  max 

If  electrons  are  introduced  into  the  wiggler  from  a  distribution  which  has 
some  energy  spread  and  is  random  in  optical  phase  we  would  expect  the  phase 
space  to  be  uniformly  occupied  as  in  Eq.  (7.28).  Such  a  distribution  would 
also  tend  to  arise  from  a  slow  diffusive  process  in  which  electrons  were 
boiled  off  from  the  top  of  the  well.  One  might  also  expect  this  distribu¬ 
tion  to  be  relatively  favorable  since  the  derivative  of  f0  lar8e  only 
near  the  top  of  the  well  where  the  coupling  to  the  radiation  is  weak. 


Since  the  Integrand  of  Eq.  (7.24)  is  non-vanishing  only  for 
J  Jffiax  we  may  take  the  orbit  Integrals  in  Eq.  (7.24)  to  be  those  for  the 
particle  at  the  top  of  the  well  i.e.  in  Eq.  (7.5)  we  should  put 

H  »  H  -  cos  ♦  -  (it  -  4>  )  sin  • 

Mr  r  r 

It  is  useful  to  take  advantage  of  periodicity  to  improve  convergence  of  the 
orbit  integral  for  particles  near  the  top  of  the  well.  We  may  thus  put 


r 


-  Q  j 


it  . 
e  +  e 


-i'l'^j  e-imfiJ^d^' [2  +  cos'l''  +  <|»'  sin^)] 

J2(HM  +  cos  ^  ^  sin  ♦  ) 


1 

r 


-1/2 

d* 


-  «K(X) 


where 


K(X) 


-iX^d'p' [2  +  cos'i''  +  V  sintr)J  *  ^ 
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+  cos  ^  sin  Yr) 


(7.29) 

is  independent  of  J  .  Note  that  the  subtraction  we  have  performed  makes 
the  Integrals  smoothly  convergent.  Equation  (7.24)  then  becomes 
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dJl^j  6(X  -  aQ)  £12|k(X)  |2/4sin  <|»r 
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But  »  ■  -  — —  and  the  sun  over  m  may  be  converted  to  an 

, 

Integral  X  d(-^)  .  Hence  we  obtain  the  growth  rate 


G(X)  - 


X|K(X) 


4  sin  tp  J 

r  max 


(7.30) 


We  have  evaluated  Eq.  (7.30)  numerically  and  find  that  the  growth 
rates  Indeed  have  the  properties  we  have  described.  The  results  are 
Illustrated  in  Fig.  7.1  for  various  .  Thus  peak  gains  of  about  G  ■  2 
occur,  typically  for  values  of  X  ~  ,4  -  .5  •  The  optimum  value  of  i|<  is 
seen  to  be  around  30°.  A  significant  feature  of  the  results  is  that  there 
is  a  range  X  <  .1  for  which  G  <  1  even  for  a  distribution  In  which  the 
well  is  filled  to  the  top.  This  implies  that  the  growing  sideband  is 
separated  by  a  finite  frequency  interval  from  the  signal  so  that  dispersive 
external  optics  may  be  applied  to  correct  the  modest  parasitic  excess 
growth. 


From  Eq.  (7.16)  we  note  that  the  unstable  region  is  given  by 
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Here  a„  ~  1  and  as  a  typical  value  is  a_  -  1(T*  ,  the  unstable  sideband 

W  8 

extends  in  the  range  10-^  >  —  >  10-^  . 
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FIGURE  7.1 

■(PARED  TO  SIGNAL.  G(  X)  IS  DEFINED  BY  EQ. 
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7.3  Discussion  and  Conclusions 

In  Section  7.2  we  have  shown  that  an  idealized  constant  wlggler 
oscillator  is  linearly  unstable  against  sideband  generation.  In  this 
section  we  will  describe  briefly  the  results  of  some  particle  simulations 
and  some  nonlinear  estimates. 

In  principle  it  appears  attractive  to  do  a  1-dimensional  particle 
simulation,  calculating  particle  orbits  and  interactions  with  the  electro¬ 
magnetic  waves.  He  have  in  fact  done  some  simulations  of  the  pendulum 
equations  (2.25)  including  the  linearized  equation  for  dtp  resulting  from 
the  presence  of  the  sideband.  By  integrating  Eq.  (7. 24)  by  parts  we  see 
that  the  rate  of  energy  transfer  by  a  particle  to  the  sideband  field  should 
be  proportional  to 

dU  ~  e2  -33  {  6(X  •  *1(++  ■♦\u|2j  .  (7.38) 


Here  E  Is  in  units  of  the  bucket  height  y  JaT  and  all  other  quantities 
are  0(1)  . 


This  result  would  evidently  be  a  bit  difficult  to  simulate!  In 
fact  the  wlgglers  which  were  simulated  were  of  finite  length  and  variable 
properties.  If  we  try  to  take  the  finite  length  into  account  by  solving 
Eq.  (7.19)  over  a  finite  Interval  in  U  the  infinite  sums  leading  to 
Eq.  (7.22)  are  replaced  by  finite  sums  and  we  find  that  Eq.  (7.38)  is 
approximately  modified  (averaging  over  the  initial  phates  of  the  particles 
or  the  sideband  waves)  by  replacing 
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V 


6(A  -  mfl)  -*•!  8l,n(X  -  •  (7.39) 

Variable  wlggler  profiles  lead  to  X  depending  on  U  as  Is  seen 
from  Eq.  (7.16)  and  It  should  be  so  interpreted  In  the  argument  of  the  6 
function  In  Eq.  (7.38).  We  have  not  found  a  simple  way  to  represent  both 
effects  simultaneously. 

The  results  of  the  simulations  so  interpreted  were  in  agreement 
with  Eq.  (7.38),  and  unstable  over-all  relative  gains  G  ■  3-4  ,  consistent 
with  the  distributions  studied,  were  observed.  There  was  however  one 
striking  result  which  has  not  been  quantitatively  analyzed.  Particles  near 
the  top  of  the  well  sometimes  gained  or  lost  large  amounts  of  energy.  This 
came  about  because  such  particles,  as  they  are  about  to  detrap,  are  In  an 
unstable  equilibrium  perched  at  the  top  of  the  well,  and  strongly  perturbed 
by  the  sideband.  For  such  particles  the  assumption  of  many  oscillation 
periods  is  clearly  incorrect.  The  fluctuations  of  energy  were  large  and  of 
either  sign  and  we  have  not  determined  whether  such  particles  produce  any 
net  effect,  either  stabilizing  or  destabilizing.  Very  good  statistics 
would  be  required  to  determine  this,  but  on  the  basis  of  our  limited 
simulations  It  seems  unlikely  that  this  effect  is  strong  enough  to  modify 
our  overall  picture. 

We  pass  now  to  a  brief  discussion  of  the  non-llnear  behavior  of 
the  excited  wave.  We  first  consider  the  situation  where  only  a  single 
unstable  mode  Is  excited  In  addition  to  the  signal.  In  this  case  we  expect 
the  wave  to  be  stabilized  by  eventual  detuning  of  resonant  particles  from 
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the  wave  as  their  energy  is  changed.  We  will  not  present  here  the  lengthy 
calculation  of  this  effect  but  will  give  a  qualitative  estimate  based  on 
Eqs.  (7*38  and  7.39). 

We  imagine  a  sideband  wave  with  amplitude  c  Interacting  with  an 
electron  beam  weak  enough  that  the  amplitude  e  is  unchanged  during  the 
interaction.  After  a  certain  distance  U  ,  which  we  wish  to  determine,  the 
resonant  electrons  will  become  detuned  and  further  transfer  of  energy  to 
the  wave  will  not  occur.  Integrating  Eq.  (7.38)  with  respect  to  U  from 
0  to  U  we  find  for  a  nearly  resonant  particle,  i.e.  a  particle  for  which 
(A  -  mJJ)U  -  0(1)  ,  an  energy  change  at  position  U 


ae 


( A-mfl) 


o 


(7.40) 


This  may  be  estimated  by  expanding  the  sine,  and  putting  (A-mft)U  ~  1  , 
while  recalling  that  in  our  units  all  parameters  are  of  order  unity.  Thus 


However,  detuning  will  occur  for  AEU  *'  1  ,  i*e.  the  resonant 

-1  /2 

electrons  will  be  detuned  for  U  >  e  .  The  width  of  the  resonance, 

i.e*  the  number  of  electrons  Involved  at  this  point  will  be  of  order 
-11/2 

U  *'  e  .  Integrating  Eq.  (7.38)  over  J  we  see  that  the  peak 

energy  transfer  rate  AE'  c  and  the  maximum  energy  transferred  to  the 
wave  in  an  infinite  length  wlggler  is  then: 
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(7.41) 
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This  contrasts  with  the  signal  wave  growth  which  continues  for  the 

Ay 

length  of  the  wiggler.  For  good  extraction,  y1  ~  0(1)  ,  this  length  is 
—1  /2 

AU  -  a  '  as  noted  earlier.  Hence  the  relative  gain  may  be  estimated  by 

/  a  vl/2 

G  *  -i--  AE  a  ^  ^  “  (—  )  .  Another  way  to  see  this  is  to  note  that  over 

£2  sat  s  \  e  / 

the  distance  in  which  it  is  amplified,  ,  the  sideband  growt':  rate  is 

only  slightly  larger  than  that  of  the  signal,  while  the  latter  continues  to 
grow  for  the  full  length  aa“l/2  #  Thus  a  single  wave  will  saturate  at  a 
very  low  level  e  ~  ag  ~  10~^  ,  since  above  this  level  its  growth  is  less 
than  that  of  the  signal  for  each  pass  of  the  wiggler. 


Unfortunately  there  may  be  of  order  10^  optical  wavelengths 
between  mirrors,  hence  10^  possible  unstable  modes  (even  neglecting  off- 
angle  propagation).  Under  these  circumstances  even  the  relatively  narrow 
resonances  for  the  saturated  single  modes  overlap  and  we  can  expect  that 
the  resonant  detuning  will  not  occur,  since  particles  will  simply  move  from 
one  wave  to  the  next  as  their  energy  changes.  In  these  circumstances  it 
seems  appropriate  to  assume  random  phases  for  the  waves,  leading  to  a 
quasl-llnear  diffusion  of  the  particles  In  the  well.  The  quasi-linear 
diffusion  equations  may  be  obtained  by  carrying  Eq.  (7.19)  to  second 
order.  Omitting  the  details  of  the  derivation  we  find: 


^“2  I  "Ij  e2  (nJJ) 
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Here  we  have  normalized  e^(X)dX  as  the  energy  in  sideband  waves 
between  X  and  X  +  dX  relative  to  the  signal.  Equation  (7.42)  is  con¬ 
sistent  with  our  physical  picture  of  the  instability,  describing  a  Brownian 
motion  of  the  electrons  in  the  potential  well  as  the  Raman  scattering  pro¬ 
ceeds.  If  we  multiply  Eq.  (7.42)  by  J  and  integrate  over  a  small 
interval  we  see  that  the  rate  of  change  in  the  number  of  action  "quanta"  in 
that  interval  is  proportional  to  m  times  the  rate  of  energy  gain  of  the 
waves  in  that  interval  as  given  by  Eq.  (7.24).  [Recall  that  f0  dj  is  the 

number  of  electrons.]  The  form  of  Eq.  (7.42)  can  be  simplified  if  we 

2  2 

assume  a  flat  sideband  spectrum,  I.e.  put  e  (X)  H  e  for  all  X  .  The  sum 
on  m  can  then  be  done  as  in  the  derivation  of  Eq.  (7.26)  to  give 
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This  is  a  diffusion  equation  with  diffusion  coefficient  of  order 

2 

e  .  Hence  a  characteristic  diffusion  distance  for  detrapping  to  occur  is 
U  ~  e  .  It  follows  that  when  the  energy  in  the  sidebands  is  greater 
than  asl/2  times  the  signal  energy  the  buckets  will  be  destroyed  before 
the  wiggler  has  been  traversed.  We  have  of  course  assumed  random  phases  In 
the  sideband  waves,  but  any  kind  of  phase  locking  which  could  restore 
stable  buckets  seems  Improbable. 


In  summary  we  find  that  the  scheme  of  decelerating  electrons 
trapped  in  potential  wells  is  subject  to  a  Raman  sideband  instability  where 
the  gain  of  the  unstable  wave  is  several  times  the  signal  gain.  It  seems 
probable  that  any  mode  of  FEL  operation  which  tries  to  suppress  energy 
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spread,  e.g.,  adiabatic  trapping,  gain  expanders,  etc.  would  be  subject  to 
similar  instabilities,  which  are  inherent  to  operation  away  from  the  peak 
of  the  linear  gain  curve. 

We  note  finally  that  rather  modest  dispersive  optics  would  be 
required  to  stabilize  a  high-Q-trapped-electron-oscillator.  If  the  cavity 
Q  were  50  it  would  have  to  be  reduced  to  20  for  the  sidebands.  This 
could  presumably  be  accomplished  by  incorporating  a  grazing  incident 
diffraction  grating  in  the  optical  system  or  by  providing  an  etalon  grating 
in  the  outer  part  of  the  mirrors.  It  would  also  be  possible  to  introduce  a 
gas  or  liquid  with  enough  dispersion  that  the  sideband  pulses  would  arrive 
back  out  of  phase  with  the  electron  micropulses,  if  breakdown  could  be 
avoided  at  high  power. 


1 1 A 


8.0 


CONCLUSION 


In  the  preceding  sections  we  have  described  and  analyzed  a  number 
of  techniques  for  making  use  of  the  properties  of  variable  parameter  wig- 
glers  to  substantially  improve  the  efficiency  of  energy  extraction  from  the 
electrons.  These  techniques  have  the  common  feature  of  requiring  the 
existence  of  substantial  ponderomotive  potential  wells  or  "buckets"  which 
are  then  manipulated  so  as  to  produce  desirable  properties.  Because  of 
this  feature,  they  all  require  high  optical  peak  power,  either  circulating 
in  a  resonator  or  emitted  from  an  amplifier,  and  are  thus  intrinsically 
high  peak  power  devices. 

A  number  of  potential  problems  have  been  identified  for  "filled" 
bucket  devices.  One  of  these  especially  applicable  to  oscillators  is  the 
potential  instability  to  parasitic  oscillations  at  frequencies  displaced 
from  the  main  frequency  by  an  amount  of  the  order  of  the  frequency  of 
particle  oscillation  in  the  trapping  buckets.  It  appears  that  some 
frequency  discrimination  device  such  as  a  diffraction  grating  may  be 
required  to  suppress  these  oscillations.  The  second,  especially  applicable 
to  amplifiers,  has  to  do  with  the  effects  of  transverse  variation.  The 
amplification  process  may  introduce  growing  inhomogeneity  in  the  optical 
field  arising  both  from  inhomogeneity  in  the  growth  rate  and  transverse 
Inhomogeneity  in  the  Induced  non  linear  index  of  refraction.  The 
assessment  of  this  problem  requires  further  work  and  clarification  before 
one  can  Judge  its  magnitude.  More  generally,  the  one  dimensional  character 
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of  our  analysis  is  a  major  deficiency.  Serious  attempts  to  design 
practical  devices  should  be  based  upon  an  analysis  which  includes  three 
dimensional  effects  more  adequately. 

Despite  the  above  comments,  it  is  our  overall  conclusion  that  the 
use  of  a  variable  parameter  wiggler  in  a  free  electron  laser  to  produce 
high  power  optical  radiation  at  selectable  frequency  and  reasonable 
efficiency  is  a  quite  promising  prospect.  The  exact  scheme  that  is  the 
best  one  to  use  (or  combination  of  various  schemes)  will  depend  upon  the 
type  of  electron  beam  that  is  available  as  well  as  the  type  of  use  for  the 
radiation.  At  the  present  time  there  are  several  groups  that  are  pursuing 
the  variable  wiggler  approach,  and  we  have  enjoyed  greatly  our  interaction 
with  them.  It  is  a  pleasure  to  thank  C.  Brau  and  R.  Cooper  of  LASL;  D. 
Prosnitz,  and  A.  Szoke  of  LLL:  S.  Mani  and  J.  Reilly  of  Schafer  Associates 
R.  Center  and  J.  Slater  of  MSNW;  and  P.  Sprangle  of  NRL  for  sharing  their 
preliminary  results  on  the  variable  parameter  wiggler  with  us.  We 
especially  wish  to  acknowledge  the  many  insights  gained  from  our 
collaborators,  V.  K.  Neil  and  R.  Novick,  without  whose  help  this  paper 
could  never  have  come  into  being. 


APPENDIX  A 


THE  RELATION  BETWEEN  ENERCY  SPREAD  AND  ENERGY 
TRANSFER  IN  THE  SMALL  SIGNAL  LIMIT 


John  Madey^  has  proved  an  important  theorem  which,  in  the  small 
signal  limit,  relates  the  phase  averaged  energy  spread  to  the  phase 
averaged  energy  change  experienced  by  an  electron  as  it  passes  through  the 
laser  interaction  region.  The  theorem  states  that 

"fY'l  (A'1) 

where  <  >  denotes  an  average  over  entry  phase  and  Y  ,  Y^  denote 

initial  and  final  values  respectively.  Y^  is  taken  to  be  fixed  so  that 

it  is  Y^  which  is  dependent  upon  the  entry  phase  before  the  average  is 

carried  out.  The  left  hand  side  is  obviously  the  phase  averaged  energy 

change,  $E  (measured,  of  course,  in  units  of  mc^).  The  phase  averaged 

L 

energy  spread,  AEg  ,  is  given  by 

AEs  e  ^f5,  -  <Yf>2  “  <(Yf  "  V2>  "  <(Yf  ■  v>2  * 

2  2 

We  shall  see  that  <(Y^  ~  Y^)>  ”  Y^)  >  vanishes  in  the  small  signal 

limit  (indeed,  it  is  implied  by  Eq.  (A-I))  so  that  (A-l)  may  also  be 


written 
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Both  Eqs.  (A-l)  and  (A-2)  are  valid  only  in  the  small  signal  limit. 


The  importance  of  the  theorem  in  the  discussion  of  the  energy 
spread  problem  in  storage  rings  has  been  indicated  in  section  3.  The  proof 
given  by  Madey  contains  a  number  of  unnecessary  approximations  and  restric¬ 
tions  which  complicate  the  exposition  and  limit  the  generality  of  the 
result.  It  therefore  seems  worthwhile  to  present  our  modified  version  of 
his  proof  here.  We  shall  see  that  the  derivation  of  (A-l)  is  not  signifi¬ 
cantly  more  complicated  than  the  evaluation  <Y^  -  Y^>  in  specific  cases, 

so  that  one  can  take  advantage  of  Eq.  (A-l)  to  obtain  <Yf  -  Y±>  from  the 

2 

much  simpler  evaluation  of  <(Y^  -  Y^)  >  . 

For  simplicity  we  make  the  (almost  certainly)  unnecessary 
assumption  that  the  transverse  canonical  momentum  is  both  conserved  and 
vanishing.  As  described  in  Section  2,  the  motion  of  an  electron  through 
the  wiggler  magnet  is,  under  these  conditions,  described  by  the  Hamiltonian 

H  -  mc2Y  ■  me2  [u2(z)  +  (  “  2aw(z)ag (z)cos  'i'  (z,t)J^2  (A-3) 

gA 

where  we  have  introduced  the  dimensionless  vector  potentials  a  ■  — -  and 

me 

V2(z)  -  1  +  a2 (z)  +  a2 (z)  (A-4) 
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aw  and  ag  are  the  vector  potential  amplitudes  of  circularly  polarized 
wiggler  and  optical  fields  respectively.  We  shall  assume  that  a  particle 
enters  at  z  ■  0  ,  t  *  0  and  exits  at  z  ■  L  .  Hie  phase  $  then 
represents  the  random  phase  of  the  optical  field  at  which  the  particle 
enters.  From  Hamilton's  equations  one  easily  shows  that 


and 


dY 

dz 


0)  a  a 
sws 


sin  ^ 


+  2a  a  cos  i|> 
w  s 


B 


cos  ♦ 


(A-6) 


(A-7) 


For  a  more  detailed  discussion  of  the  equations  of  motion,  the  reader  is 
referred  to  Section  2. 


We  are  now  in  a  position  to  undertake  the  derivation  of 
Eq.  (A— 1 ) .  First  write 


y  -  \  +  Yj  +  y2+ 


(A-8) 


where  is  first  order  in  as ,  Yj  second  order,  etc. 
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From  (A-6)  ve  have 
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a  (z')a  (z')sin  ^  (z' ,t(z' )  )dz' 
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J  Y2(z')  -  W2(z')  +  2aw(z')  as(z')  cos  *  (z'.t(z')) 

(A-9) 


where 
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-  t  +  t,  +  (A-l 1) 

o  1  •  •  • 

and  t  ,  ti,  etc*  also  refer  to  zero,  first,  etc.  order  in  a„  .  we  have 
o  1  “ 

immediately 
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where  =  ^(z'  ,tQ(z')  J,  \?Q  ■  1  +  ,  and  in  general,  the  primes  indicate 

that  the  argument  is  z'  (and  double  primes,  that  the  argument  is  z"  ). 
Using 


<sin  4»'>  -  0 
o 


(A- 13) 
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and 


<sin  <l>'  sin  ty">  -  -j:  cos  (<i>'  -  V") 
o  o  2  oo 


(A- 14) 


yields 


<Yj>  «  0 


(A- 15) 


and,  to  lowest  nonvanishing  order  (second  in  ag  ) ,  using 

<(Y,  -  Y,)2>  -  <  Y2  >  _  ,  also  yields 

1  t  1  2  “  L 


2  L  L 


<(V  V  > 


(u  r  r  a  a 

t>  /“  tH- 


a"a" 

W  8 


O  O 


Y2  -  p'2  /Y2  -  p"2 

Ti  o  r  i  o 


COS(<il'-  <J|"  ) 

o  o 


(A- 16) 


Expanding  Eq.  (A-9)  one  obtains 


'  C  J  *i 


I  2  ,2 

Y,  -  P 
i  o 


Y!  sin  V  - 
1  o 


u>  a'a' 
s  w  s 


Y2  -  P'2 

i  o 


t'  cos  dz' 

1  o 


»  zr 

*1 


ar  a*  ^ 

a  a 

W  8 


— - - sin  cos  V  dz' 

(r  -W2)3/2 

O  10 


(A-17) 


We  write  Yj  for  the  first  line  of  Eq.  (A-17)  and  note  that  since 


<sin  V  cos  <!>'>  -  0,  <Y,>  -  <Y„>  . 

o  o  2  2 
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Next,  expanding  Eq.  (A-10),  and  using  Eq.  (A-7)  yields 


and  writing  t^  -  tJa  +  t^ 


(A-18) 


(A-19) 


t 


la 


dz" 


(A-20) 


In  writing  Eq.  (A-21),  we  have,  in  anticipation  of  what  is  to  follow,  made 
use  of  the  identity 


Substituting  Eq.  (A-12)  in  Eq.  (A-20)  and  interchanging  the  order  of 
integration  we  obtain 
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dz" 


(A-23 


Now  substitute  Eqs.  (A-12),  (A-21)  and  (A-23)  into  Eq.  (A-17)  to  obtain 


(A-2- 


Then  phase  averaging  and  supplementing  Eq.  (A-14)  with 

<cos  cos  cos  (♦*  “  W)  (A-2! 

o  o  z  oo 


<cos  sin  »!>">  - - ^sinfiK  -  ♦'') 


(A-2( 


yields,  to  lowest  nonvanishing  order  (second  In  aa) 

s 

<Yf  -  V  -  "V  z-L 


dz‘ 


cos(<K 

o 


) 


-4  <<Tf  -  v2> 


(A-27) 

(A-28) 


In  going  from  Eq.  (A-27)  to  Eq.  (A-28)  we  have  made  use  of  Eq. 
(A— 1 6)  and  of  the  symmetry  in  (z',z")  of  the  integrand  of  Eq.  (A-27). 
Eq.  (A-28)  and  Eq.  (A-l)  are  the  same  so  the  proof  Is  now  complete. 


It  should  be  added  that  an  identical  theorem  holds  for  the  Yariv, 
Shih^2  device.  The  proof  is  similar  but  simpler. 

Madey  has  also  derived  Eq.  (A-l)  by  means  of  a  quantum  mechanical 
argument  in  which  stimulated  and  spontaneous  emissions  are  related  to  one 
another  and  to  Eq.  (A-l).  In  a  similar  spirit,  we  offer  an  even  simpler 
proof  based  upon  the  principle  of  detailed  balance.  Let  P^.Yj)  be  the 
probability  density  In  that  an  electron  which  enters  at  y  *  Y^ 

leaves  with  Y  “  Y2  •  PrlnclP^e  detailed  balance  states  that 

P(Y1,Y2)  -  P(Y2,Y1)  .  (A-29) 
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Assuming  that  the  wiggler  magnet  with  optical  field  satisfies 
Eq.  (A-29) ,  we  have,  with  “  Y2  ”  ^1  * 

POfj.Yj  +  ay)  -  p(y1  +  ay,y2  +  ay  -  ay) 

-  Pttj.Yj  -  Ay)  +  ^  P  (Y1,Y1  -  Ay)  Ay  (A-30) 

Eq.  (A-30)  holding  only  in  the  small  signal  limit.  Now  multiplying  Eq. 
(A-30)  by  Ay,  integrating  over  Ay,  and  using 

< (Yf  -  Y±)>  -  /P[Y1.Y1  +  AY) AY  d(AY)  »  -  J?{ Yj.Yj  -  Ay)Ay  d (Ay) 


and 


m 


ay 


(Ay)  d ( Ay) 


a 

■3y7 


7  /' 


p(VYt 


-  ay)  (Ay)  d (Ay) 


) 


2 

> 


we  obtain  Eq.  (A-l)  again. 

It  is  perhaps  worth  emphasizing  that  Eq.  (A-l)  demonstrates  that 
for  finite  wigglers  in  the  small  signal  limit,  zero  energy  spread  implies 
zero  energy  transfer.  To  see  this  we  note  that  Eq.  (A-16)  may  be  written 
in  the  form 

u)2 

<(Yf  -  Yt)2>  F*F  (A-31) 
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with 


(A-32) 


The  vanishing  of  <(Y^  -  Y1>^>  therefore  Implies  that  F  vanishes.  For 

finite  L,  cannot  be  singular  so  ^  <(y  -  Y  )^>  must  vanish  as 

dY^  oij  f  1 

well. 


As  a  simple  example  of  the  application  of  (A-l)  we  consider  the 
case  of  a  variable  parameter  wlggler.  Equation  (A-5)  specialized  to  zeroth 
order  yields 


di|> 

O  | 
ST  m  kw 


2Y2c 


(A-33) 


and  for  simplicity  we  consider  the  case  studied  by  Brau^S 


k  -  k  +  k'z  (A-34) 

w  wo  w 

with  a  ,a  ,u  constant  and  v/y^_  y2  ;  «  .  Because  of  Eq.  (A-31)  we  may 

w  s  1  o  1 

choose  the  Integration  constant  arbitrarily  in  determining  ♦  from  Eq . 
(A-33).  Choosing  this  constant  conveniently,  we  obtain 


<(Yf  -  Yt)  >  - 


w^a^a^ 

s  w  a  F  F* 

2c2Y?  11 


(A-36) 
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where 


■I 


dz  exp  1 


k' 

w 

T 


(z  - 


zo>‘ 


with 


(A-37) 


z 

o 


(A-38) 


The  point  zq  corresponds  to  d^/dz  as  given  by  Eq .  (A-33)  .  It  is,  of 

course,  the  phase  matching  point  for  given  w  ,  Y  when  It  occurs  within 

s  1 

the  wiggler.  It  Is  apparent  from  Eq.  (A-37)  that  the  dependence  on  zQ  is 
identical  to  that  which  occurs  for  the  single  slit  Fresnel  intensity 
problem,  where  zQ  Is  the  displacement  of  the  observation  point  from  and 
parallel  to  the  slit  edge  and  1^/k^  its  distance  from  the  slit  plane. 

The  constant  wiggler  corresponds  to  the  Fraunnofer  limit  and  hence  to  the 
ca8e  *;l2  <  1  exPress  Fj  in  terms  of  standard  Fresnel  integrals  we 

define 


F2<y) 


2 

z 


dz 


C  [y  /k'/r  )  +  IS  [y  /k'/v  ) 

W  w 


(A-39) 


Then 

F i  -  F2  (zq)  +  F2  (L  -  zo) 


(A- 40) 


end  Eq.  (A-l)  yields  immediately 


2  2  2 
w  a  a  dz 

S  W  8  O  _ 

<yf  -  Y  >  -  - -7T-  Re 

f  °  ,  2c2Y2  dYi 


3  2  2  2 
to  a  a  n 
s  w  s 


2c3Y,V 

i  w 


dTFi 
o  / 


.  kw,  .2  .  kw  2 

-i  —  (L  -  «  )  -i  —  *0 

e  -  e 


(A-41) 


C^_ (*  )  +  F  (L  -  z  )} 
2  o  2  o 


Equation  (A-41)  reveals  the  energy  transfer  to  be  a  rapidly 
oscillating  function  of  zQ  .  While  its  derivation  is  simple, the 
discussion  of  its  behavior  is  complicated.  For  z  well  within  the 


wiggler,  F2(zq)  +  F2(L  "  Zo) 


—  and  hence 


Re  -j— ^  F, 
dz  1 
o 


I - /  k'  k'  ' 

/_n_  f  w  ,T  n21,'i  r  w  -  2  fT  y 

lw  C0S  1  r  (L  ■  Z0)  -  4  1  -  cos  (  T  Zo  -  4  ] 

w 


^|Tsin  k^  sin  (  -f  (ij  +  zj)  ) 


where  Zj  »  zQ  -  L/2  is  the  distance  of  zQ  from  the  wiggler  center.  The 
energy  transfer  is  seen  to  oscillate  symmetrically  about  zero  with 
decreasing  period  and  non  decreasing  amplitude  as  zQ  moves  towards  the 
center  of  the  wiggler.  An  energy  spread  in  the  incident  beam  induces  an 


I* 


128 


For  typical  values  of 


average  in  z 


2^ 

over  a  range  &zq  jpr-  — — 

w  ’i 


6^ 


and  the  other  parameters  the  average  extends  over  several  periods 
(~  5)  near  the  center.  For  such  a  beam  the  oscillations  will  decrease 
strongly  in  amplitude  as  one  moves  towards  the  center. 

The  situation  near  one  side  of  the  wiggler  is  different.  For  z0 
negative  or  “  zero  we  write 


dz 


r  k<  o 

.  w  2 

/  \  ""^2  2o 

It' 

,  w  (  2  ,T  « 2  \ 

V^o  "  (L  '  Zo)  J 

F  -  - 

IF  (z  )  +F„(L  -  z  )]  e 

1  -  e  2  °  ° 

1 

\  2  o  2  o  ) 

- 

The  first  factor  in  square  brackets  is  non  oscillatory  in  the  range 
considered.  The  second  factor  is  even  more  rapidly  oscillating  than  before 
but  with  average  value  equal  to  1.  Again  an  energy  spread  will  strongly 
damp  the  oscillations.  We  therefore  omit  the  oscillating  term  to  obtain 


Re 


7T  rl  --('/*^f)s  {*oJ *V’) 


where 


g(x)  = 

j  -  C(x) 

*  2 

cos  x  + 

\  -  S(x) 

■  . 

1 

2  +  4. 142x  +  3.492x2  +  6.670x3 


,  *  2 
sin  -  x 


for  x  >  0  . 
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For  x  <  0  ,  g  reaches  a  peak  of  '  1.3  at  x  ""  -.74  and  then  commences  to 
oscillate  at  non  decreasing  amplitude  and  decreasing  period,  being 

sin  (^x2  +  ?-)  for  |x|  large.  The  situation  at  z  >  l  or  “  L  Is 

L  4  O 

similar,  but  with  the  sign  of  the  effect  reversed.  A  qualitative  picture 
of  the  behavior  as  a  function  of  <Y^>  that  results  after  energy  averaging 
is  shown  in  Fig.  A-l. 
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APPENDIX  B 


EVALUATION  OF  PHASE  AREA  DISPLACEMENT  ENERGY  TRANSFER  AND  ENERGY  SPREAD 


We  may  obtain  a  formal  result  for  Ay  by  rewriting  Eq.  (6.13) 


^f  7*  ( 

Ay  *  -  ^ ^  */ a"  sin  ^  J  J  dX  <  e 

o  I 


(cos  ^ 


l>  +  ij’sin 


-X^Jh  +  (cos  ^  sin 


(B-l) 


where  here  H  ■  -cos  i{i  -  ^sin  «  and  we  focus  our  attention  upon  the 
turning  point  which  lies  between  the  minimum  at  and  maximum  at 

as  shown  in  Fig.  6«2.  The  range  of  within  these  limits  is  given  by 
i{>  <  ^  <  ir-^r  ,  where  ^  satisfies 


cos  +  ty.sinty  -  cos  (tt  +  ^)-(7T4-^)  sin  ^ 

1  1  r  r  r  r 


r 


and  is  given  by  ^  -  (ir  -  </4*|  ^1  -  ^r)  for  small  |  .  Since 

is  negative  here,  we  find  it  convenient  to  write  (|>  ■  -|v  |and  replace 
by  ^  -  n  ,  so  that  -  ^/4 n |  ^  |  <  ^  -  1^1  <  o  .  Making  these 

changes  and  also  shifting  the  variable  of  integration  by  n  we  obtain 


Ay  =  2  /2A/it  sin |  \|»  | 


dX’exp  X  [cos  cos  ^  +  (’*'  -  4^)  sin|^|] 


-exp 


[X2(*  - 


sin  |  V 


(B-2) 


Next  we  write 


exp  X  (cos'*'  +  sinl^l)  d^ 


±J 

n*=0 


C  -  2"n 


exp  X  (cos  'J'  +  sin|'i'  |  ) d v 


-  2*(n  + 


1) 


d4> 

\  -2  it 


exp  X^(cos  'j'  +  4>  sin  |  4^  |  )  [exp(-2TX^sin|  4^1 )  ]n 


n“0 


2 

exp( xX  sin | 4*  | ) 

2 

2  sinh(*X  sin|  4»  | ) 


2 

exp  X  (cos  4*  +  41  sinl^l)  d'l' 
2w 


(B-3) 
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and 


» 

I  exp  X2 (cos  *  +  ♦  sinH»rl)d*  *  /  exp  X2(cos  ♦  +  *  sinl^Dd^ 

^’-2*  -v 

+  exp(-2nX2sin|  i>  I)  /  exp  X2(cos  Y  +  Y  slnl^l)  dV 

r  4* 


(B-4) 


Substituting  (B-4)  into  (B-3)  and  (B-3)  into  (B-2)  we  obtain 


&Y  *  2  \/2  A/X  sin  |  ♦ 


r  -* 

1  |  I  dX  e 

r  i 


2  .  ' 
cos  v. 


2  sinh  [xX  sin  1^1  ) 


-nX2sin|*rl  -X2sinKrl  -  *)  X2cos* 

•{,  I 


nX2sin|  V  |  f^l  -X2sin|  V  |  X2cos  » 

+  e  r  /  e  r  f  e  d* 

iv 


X  cos 


X  sin | I 


(B-5) 


which  we  rewrite  as 
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Ay  -  2  / 2A/TT  sin|  |  /  dX 


■'/ 


,2 

-X  cos 


r2sinh  TtXZ  (s In  |^»r  |  ) 


cosh  (ttXZ 


sin| trl ) 


X  cos 
e 


* 

d* 


+  sinh(nX  sini^l) 


■ 

-u 


d^  e 


-Xzain|  ^K^- 


.2 

X  cos 

f  j 


X  sin |  i|»r  | 
(B-6) 


This  Is  of  course  a  bit  complicated.  However  we  may  perform  the 
Integral  over  \|>  In  the  limit  <J/'2,|<J<r|  +  0  to  yield,  after  some 
manipulation, 


+ 


e 


x2/2 


(B-7) 


where  B  =  Xj ^  -  |i|i  ||  and  I0  Is  the  usual  Bessel  function  of  Imaginary 
argument.  This  Is  still  too  hard  to  do  exactly.  Splitting  it  Into  ranges 
we  have  for  the  Integrand  if 
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>  A 


A2U 


I  (X2)  e"X 

o 


1  + 


(*J-I*rl)2|^"  coth  .X2|*rl  -/?>  x 


*"d  11  1  ’ 


i  (x2,  ,-*2  .T<*r'*r"‘ 

o 


A2  2 

i*j> 


2  it  e 


-2«^l*rl  yi7e  2  VTf  '  r 


A  |  V£'-l  ^  II 


1 


A2Krl 

(B-8) 


The  integrations  may  be  done  approximately  to  yield 


Ay  m  2  yJlAh 


00 


[i  (A2)  e"X  -  lj 


dA  -  °  ■' - = - -  + - ~  (*f~  l^rl>2  *n 

A2  2  y/T " 


v/sin|  ^  | 


(B-9) 


with  El(z)  the  exponential  integral.  We  observe  that 


/ 


-A  2 

e  I  (A2)  -  1  - - 

d  A - 1 - -  -2  / 2/ f 
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mm 


Note  that  the  result  diverges  unless  -  \/ 4 it  1 |  <  -  Ill'll  <  0 

which  is  the  allowable  range  for  turning.  The  derivation  is  correct  only 
in  the  limit  of  very  small  |i|>  |  •  The  logarithmic  terms  are  exact  but 
errors  of  order  ( in )~*  are  to  be  expected.  We  have  somewhat  arbi¬ 
trarily  chosen  them  so  that  the  leading  term  goes  to  zero  for  sinful*  1  • 

Noting  that  the  range  of  ♦  '  -1^1  is  from  _  y ^  |  ^  |t0  0  , 

we  put  -  1^1  )  ■  -^4 tt x  •  Then  as  a  rough  approximation  we  have  with 

o  <  X  <  1  , 


ay 


4  yrl-i 


1T^r 

___ 


x2  I  An  sinN»rl|  +  |An  xl  +  An  1 1-Xl 


(B-10) 


and  doing  the  averages  over  2xdx  we  find  the  r.m.s.  spread 


AY 

rms 

Ay 


[  -  1  sin| -  3  )  +  3  I  An  sin|^r  ||  +  I  An  sinUJ 


(B-l 1) 


Note  that  AY  -  -  ^  is  exactly  what  is  obtained  in  the  0  limit 

by  taking  the  phase  space  displacement  at  the  bump  as  expected.  The  terms 
proportional  to  | An  x!  and  An ( 1 — X I  in  Eq.  (B-10)  arise  of  course  from 
particles  which  just  graze  the  top  of  the  reflecting  wells. 
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A  further  energy  spread  Is  to  be  expected  from  the  fact  that  the 
bucket  is  of  course  not  moved  up  from  -*>  to  +“  but  has  a  finite  energy 
traverse  2Ay  .  Thus  in  Eq.  (6.12)  if  we  look  at  the  end  point  contribu¬ 
tions  to  &Y  we  see  that  we  may  expect  a  spread  from  the  two  end  points 


Ay  «  cos  ij<  A  ^2  . 

-w 

Ay 

Eliminating  A  by  using  Eq.  (B-8)  we  find  for  the  r.m.s. 
contribution  to  end-loss 


Ay 

rms 

'Ey 


(B-12) 


The  spreads  (B-ll)  and  (B-12)  are  presumably  statistically 
independent  so  they  should  be  added  on  in  the  r.m.s.  sense.  For  small 
resonant  angles  we  may  approximate  the  slowly  varying  terms  in  Eq.  (B-ll) 
and  combine  it  with  Eq.  (B-12)  to  give  a  simple  formula  for  the  ratio  of 
spread  to  gain 


Ay 

rms 

Ay 


(B-13) 
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APPENDIX  C 

THE  LIMITATION  ON  ELECTRON  BEAM  RADIUS  IMPOSED  BY  EMITTANCE 

All  of  the  discussion  in  the  main  body  of  this  paper  has  been 
carried  out  under  the  assumption,  P^  *  0  .  Any  real  electron  beam  will, 

however,  have  a  distribution  of  transverse  velocities,  so  that  at  best  one 
can  only  expect  pj/pz  t0  small*  While  we  are  not  prepared  to  discuss 
the  effects  of  this  spread  completely,  there  is  one  important  effect  which 
can  be  readily  estimated  and  which  provides  a  basis  for  estimating  the 
magnitude  of  spread  which  can  be  tolerated.  This  effect  is  the  spread  in 
the  axial  velocity  vz  which  a  spread  in  implies  for  a  monoenergetic 
beam. 


We  assume  an  azimuthally  symmetric  beam.  The  phase  area  in  x  , 
•gi  which  contains  the  beam  is  equal  to  that  in  y  ,  ijZ  and  is  denoted 
by  ne  .  The  quantity  e  is  referred  to  as  the  beam  emittance.  We 
therefore  take 


P 


1 


(C-l) 


as  an  appropriate  measure  of  the  variation  in  p  (  rfi  is  the  electron 
beam  radius  as  in  Section  4).  We  take 


y 


2 
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1/2 


(C-2) 


to  be  fixed,  so  that  the  deviation  of  ?z  oscillates  with  the  wiggler 

period.  Since  this  period  is  short  compared  to  the  synchrotron  period,  it 

is  appropriate  to  average.  Hence  we  obtain  (assimlng  6  ~  1  ) 

z 


<6p  > 
z 

P" 

z 


<63  > 
z 

- B 


(C-3) 


We  require 


<66  > 
z 


6y 

max 


(C-4) 


where  ^max  is  bucket  height  defined  by  Equation  (2.56).  Combining 
these  relations  yields  the  inequality 


For  linear  accelerators  the  following  relation  has  some  theore- 

/ 

tical  and  considerable  experimental  support^ 


K__  T  1/2 

Y6 


e 


(C-6) 


where  K  is  referred  to  as  the  emittance  constant.  Ig  is  the  peak 
electron  beam  current  averaged  over  a  macropulse.  Hence  we  obtain 


r 

e 


> 


I 

e 


1/2 


K 


( C—  7 ) 
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All  of  the  factors  which  appear  in  Equation  (C-7)  can  be  related  to  the 
"standard  design"  parameters  discussed  in  Section  4.  Thus  using  Equations 
(4.25)  and  (4.26)  we  obtain 


y/2  ^3/2  (2-^)5/2  csc3/2*r 

n^1-^)  (2-2 vV>  r 


1/2 


•  K 


(C-8) 


which,  for  the  numerical  example  following  Eq.(4-26),  takes  the  form 


(C-9) 


Since  Ig  refers  to  the  peak  current  in  a  micropulse,  the  factor 

Ie/Ie  varies  from  1/25  to  1/100,  depending  upon  design.  A  typical  "good" 

value  for  K  is  .3  cm  -  (kA)”*^  •  Assuming  this  value  for  K  we  see 

that  the  choice  re  *  1.24  mm  used  in  the  numerical  examples  satisfies  C-9 

for  I  /QI  <  1/163  . 

e  e 


In  the  case  of  amplifiers,  the  peak  current  required  is  so  large 
that  it  is  probably  necessary  to  use  an  induction  accelerator  instead  of  an 
RF  llnac.  Then  the  factor  I^/l^  is  simply  unity.  In  addition  Q  is  of 
course  also  unity.  Actual  numbers  depend  upon  gain  required  and  wiggler 
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profile  design  details,  but  for  our  example  of  Figure  4.8,  we  find 
rg  >  5.6  K  cm  .  Reference  to  Equation  (4.20)  suggests  that  an  order  of 
magnitude  reduction  in  K  would  greatly  enhance  the  practicality  of 
amplifier  applications. 

It  should  be  mentioned  that  emittance  also  affects  the  beam 
transport  problem^.  Furthermore,  to  the  extent  that  one  depends  upon  the 
wiggler  to  provide  focussing,  there  is  a  connection  between  emittance  and 

the  wiggler  factor  o  (Eq.  4.10a)^»^.  These  considerations  have  not 

w 

been  taken  into  account  in  deriving  Ec .  (C-8) . 


APPENDIX  D 


NEGLECT  OF  SPACE  CHARGE  FORCE 


In  this  appendix  we  estimate  the  importance  of  space  charge  forces 
in  the  discussion  of  longitudinal  motion.  We  do  not  consider  transverse 
space  charge  effects. 


by 


The  ponderamotive  force  on  an  electron  has  a  z  derivative  given 


mu  cosy 
s  r 


(D-l) 


at  the  bottom  of  the  bucket.  K  may  be  thought  of  as  the  bucket  spring 

P 

constant.  The  negative  spring  constant  arising  from  space  charge  forces  at 
the  bottom  of  the  well  is  given  (approximately)  by 


K 


sc 


2*f.  ep 
b  e 


2fbTee 


2 

r  c 
e 


(D-2) 


where  p  is  the  electron  beam  charge  density  and  f  is  the  trapped 

e  b 

fraction.  Hence 


2fb  V  e 


2  2 

emu  a  a  r  cosy 
s  w  s  e  r 


(D-3) 


We  consider  space  charge  forces  to  be  unimportant  when  this  ratio  is  small. 
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From  2.40  and  2.41  in  the  standard  model 


V2-V 


s  2k  L(l-n.  )sin* 
w  b  r 


(D-4) 


From  4.24  and  4.11  (taking  n  »  f  n  ) 

e  b  b 


i  3 

1  me  b 


(2-nbV 


0  2  e  f  2  2 

b  (l-nb) sin  *  QY£ 


and  from  4.4  at  z  «  L  ,  where 


M>) 


l  is  largest, 


(D-5; 


(D-6) 


Combining  the  above  relations  in  the  standard  model  we  obtain. 


Ksc  .  (2~2Vnb')  <2~V 

K  "  2  3 

p  8QY  ( 1  — n,  )  sin^  cos^ 


(D-7) 


For  the  values  of  Y^  which  we  have  in  mind,  this  is  always  a  very  small 
number.  For  the  numerical  example  giver,  after  Eq.  4.26, 


K „  /  «  5.8  x  10"4  /  Q  . 

sc  p 


For  the  case  of  an  amplifier  the  Q  factor  is  of  course 
missing.  The  space  charge  effect  is  relatively  strongest  at  the  input  end 
because  the  ponderamotive  force  is  weakest  there.  The  numerical  factors 
depend  upon  design  details.  For  the  amplifier  of  Figure  4-8  we  find 
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CD-8) 


so  that  space  charge  effects  can  be  significant  for  smaller  values  of  Y^ 
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APPENDIX  E 


THE  OPTICAL  KLYSTRON 


In  order  to  improve  the  efficiency  of  the  constant  parameter  FEL 
in  storage  ring  operation  it  is  logical  to  think  of  bunching  the  electrons 
before  sending  them  through  the  wiggler.  One  could  speculate  that  by 
properly  choosing  the  phase  of  the  modulation,  energy  spread  could  be 
greatly  reduced^.  To  study  this  we  consider  a  system  of  two  wiggler 
magnets  (a  velocity  buncher  and  a  radiator)  separated  by  a  long  drift 
epace  (Fig.  E.l). 

0 

!<-<><  -  L  ><  /»  z 


Fig.  E.l  Schematic  Representative  of  Optical  Klystron 
Configuration 


This  is  a  special  soluble  case  of  a  variable  wiggler,  where 
bunching  becomes  large.  The  electron  motion  is  given  by  Eq.  (2.44) 


d± 

dz 


-A(z)«Y 


(E-l) 


d^Y 

dz 


C(z)  sin  <i>  . 


(E-2) 


149 


•I 


Here  A  and  C  are  defined  in  the  main  text.  For  our  case, 

since  A  «  L  ,  we  may  take  A  -  0  for  z  <  A  ,  and  z  >  &  +  L  ,  and 

A  “  const  in  the  drift  space.  We  take  C  ■  Cj^  in  the  first  magnet  and 

C  »  C2  in  the  second.  Consider  an  electron  with  6y  ■  6y  ,  ip  »  ^  at 

oo 

z  «  0  .  Then  at  z  «  A 


♦  -  ♦  ,  -  6Y  +  C. ASin  <j> 

o  o  1  o 


at  z  =  A  +  L 


^  “  A  J^6Yo  +  C  x  As  in  L  +  X  , 


where  we  have  allowed  for  a  phase  shift  x  in  the  second  magnet  and 


5y  ■  <$Y  +  C,  Asin  <|i 

o  1  o 


and  finally  at  z  =  L  +  2A 


<5y  -  «Y  +  C  Asin  *  +  CAsml*  +X-AL/6Y  +C,Asin 

o  1  o  2  [  o  \  o  1  o/J 

ao 

-  ^Y  +  C  Asm  *  +  C,aV  j  (AC, LA)  sin  |  ♦  +  x  -  AL«Y  -  nY  1  , 
o  1  o  2  ^  n  1  L  o  o  oj 


(E-3) 


with  Jn  being  the  nc^  order  Bessel  function. 


We  may  now  average  over  initial  phase  to  obtain 

o 
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<6y  -  6y  C0iJ,(AC,Li)  sin  (x  -  AL ) 
\  o/  2  l  l  o 


(E-4) 


<(6y  -  fiYo)2^>-  ^  C2*2  C2A2[l-cos  2(X-AL6Yo)J2(2AC1L4)] 


+  C. C„ *cos( X  -  AL6y  )  <  J  (AC, Li)  -  J„(AC,LJl) 

12  OIO  1  21 


(E-5) 


<(«Y  -  «V2> 

If  we  now  minimize  the  ratio  - 5—  with  respect  to  Co/Ci 

/6y  -  6yo>2 

holding  other  factors  fixed  we  have: 


(<5y-<5y  ) 


<(fi Y  -  6V/ 


1-cos  2(x-AL<5y  )  J_  (2AC.LE)  -  cos2(x-AL<5y  )  J  -J  2 


JjSin  (X  -  AL6Yq) 


L  i-j2<2B)-Cj0(P)-j2 ( B) )2  2J2(26)+[jo(e>-J2(e))' 

2  J2(g)  sin2(X  -  AL«Yo)  jJ(8) 


with  8  ■  AC^Lfc 


Since  1  -  J2(28)  -  jjo(6)  ~  J2(8)j2  >  0  ,  the  optimum  value  for 
the  phase  will  be  such  that  sin  (X  -  AL5Yo>  ■  1  ,  i.e.  Cj/C2  “  0  . 

Hence,  finally  for  the  optimum  energy  6y  ,  the  ratio  of  spread  to  gain  is 


given  by: 


<(»Y  -  SY0)2)  , 

<(«,  .  6Yo)>2  '  2 


/(«Y  -  «Y  )2)> 

The  ratio  Is  optimal  for  6  -  2.2  where  — - - r-  «*  2  . 

<(«Y  -  «Yo)>2 

This,  at  first  sight,  may  appear  advantageous  for  finite  B  operation. 
However  we  may  note  that  for  small  values  of  B  the  moments  given  by 
Eq.  (E-4)  obey  the  Madey  theorem: 


1  +  J2(2B) 


J*(S) 


(E-6) 


<(6Y  -  6tq)>  =^-4-  *  6V2) 

o 


1  Cl4 

while  as  B  is  Increased  the  bandwidth  ^  -~jp  also  becomes  very 

small.  For  storage  ring  operation  the  ratio,  bandwidth  x  mean  energy 
loss/ (energy  spread) optimised  with  respect  to  C2/C^  for  &Yo  at  the 
peak  gain  point  (sin  (x  -  A£6Yo>  *  -1),  Is  presumably  more  relevant.  The 
optimised  ratio  Is  found  to  be  "  J^(B)  /  B  so  that  low  B  operation  Is  In 
fact  preferable.  Similarly,  use  of  the  Fokker-Planck  equation  to  determine 
a  steady  state  solution  from  Eqs.  (A-4)  and  (A-5)  shows  that  even  for  large 
B  ,  the  steady  state  is  characterized  by  an  emitted  power  which  Is  a  small 
fraction  of  synchrotron  radiation  (and  a  decreasing  function  of  B  ). 
However,  that  is  only  an  indicative  result,  not  a  conclusive  one  since  for 
large  B  ,  the  Fokker-Planck  equation  is  not  valid  as  higher  derivatives 


ilHi  become  Important  as  the  bandwidth  becomes  small.  Pending  a  presumably 


3Y 


numerical  solution  of  the  Integral  equation  for  the  distribution  function 


In  a  storage  ring  at  large  &  ,  we  can  only  say  that  the  outlook  is  not 
promising.  In  particular  the  result  (E-6)  shows  that  it  is  not  possible, 
even  for  a  5-function  input,  to  completely  eliminate  energy  spread. 


i 

i 
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